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ABSTRACT 


A  wind  tunnel  model  was  used  to  study  stack  plume  dispersion 
at  a  surface-mining  site.  The  interaction  between  plumes  and  the 
wake  flows  behind  the  stacks  and  large  terrain  obstacles  (dikes) 
was  investigated. 

Mean  flow  Reynolds  number  was  about  25,000  times  smaller  in 
the  wind  tunnel  than  in  the  full  scale  atmospheric  flow,  affecting 
the  dynamic  similarity  of  the  model.  Stack  wake  flows  were  found  to 
model  downwash  effects  properly  despite  Reynolds  number  mismatch,  as 
indicated  by  measured  values  of  base  pressure  coefficient,  -Cp^  -  0.86, 
for  cylinders  at  model  conditions,  which  were  close  to  values  from 
Roshko  (1960)  for  cylinder  flows  at  full  scale  Reynolds  numbers. 

However,  model  dikes  were  found  to  greatly  exaggerate  the  size  of 
recirculation  zones  and  influence  of  the  dike  wake  flows  on  downwind 
velocity  and  turbulence  levels.  This  problem  was  corrected  using  a 
deflector  vane  mounted  on  the  model  dike  crests. 

A  combined-rise  formulation  was  developed  which  predicts 
both  the  momentum  and  buoyancy  effect  on  plume  rise.  Momentum  rise 
entrainment  constants,  predicted  from  a  simple  model  from  Wilson 

(unpublished)  were  found  to  yield  accurate  momentum  rise  predictions 
using  the  combined-rise  formulation.  Buoyancy  rise  entrainment  con¬ 
stants  $2  were  found  to  vary  realistically  with  changing  flow  conditions. 

The  momentum  flux  ratio  d)..  was  found  to  be  well  correlated  with 

M 

stack  downwash  effects,  providing  a  means  of  prediction  for  this 


(v) 


phenomenon . 

Dike  wakes  were  found  to  increase  measured  ground  level 
concentrations  (GLC)  from  plumes,  with  a  non-linear  degree  of  severity 
with  dike  height.  The  Gaussian  dispersion  model  was  found  to  be 
adequate  in  predicting  GLC.  A  correction  model  was  proposed  to  account 
for  dike  effects  on  GLC  through  increased  vertical  plume  spread  o 
in  the  Gaussian  model. 
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CHAPTER  I 


INTRODUCTION 

Plume  dispersion  from  tall  industrial  stacks  has  been  studied 
in  detail  in  recent  years,  largely  due  to  the  increasing  stringency 
of  pollution  control  legislation  and  the  growing  concern  for  all 
industrial  projects  to  be  accountable  for  their  impact  on  the  environ¬ 
ment.  The  designer  of  a  plant  must  be  able  to  predict  the  maximum 
levels  of  pollutant  concentrations  that  are  likely  to  occur  at  ground 
level  or  at  other  points  of  interest. 

In  this  study,  two  special  problems  encountered  in  the  pre¬ 
diction  of  pollutant  concentrations  from  stacks  are  investigated. 

Both  deal  with  the  interaction  of  wakes  behind  man-made  flow  obstacles 
with  stack  plumes.  The  first  wake  effect,  common  to  all  stacks,  is 
the  interaction  of  the  plume  gas  with  the  wake  generated  by  the  stack 
body,  acting  as  a  flow  obstacle  in  the  atmospheric  boundary  layer. 

This  small  wake  flow  can  have  a  profound  effect  on  the  dynamics  of  the 
plume  motion  because  of  its  close  proximity  to  the  plume  in  the 
critical  early  stages  of  rise.  Under  certain  conditions,  the  low 
pressure  in  the  stack  wake  can  pull  the  plume  down  from  its  normal 
rise  trajectory  and  cause  abnormally  severe  pollution  effects.  This 
phenomenon  is  called  stack  downwash  and  has  been  observed  by  many  in¬ 
vestigators  such  as  Slawson  (1978)  and  Fay,  Escudier,  and  Hoult  (1969). 
Briggs  (1969)  gives  a  history  of  the  study  of  stack  downwash,  which 
dates  back  to  the  1940's.  A  simple  criteria  from  that  era  for  avoiding 


1 


2 


downwash  was  to  ensure  that  R>1.5,  where  R  =  w  /U  is  the  ratio 

s  s 

of  the  stack  efflux  velocity  to  the  windspeed  at  stack  height.  A  more 

general  criteria  is  developed  in  this  study  which  uses  the  momentum 

2  2 

flux  ratio,  <J>.,  =  ap  w  /p  U 
M  s  s  a  s 

It  is  necessary  to  understand  how  stack  downwash  affects  plume 
rise  before  investigating  the  effects  of  other  plume  disturbances, 
because  two  or  more  plume-wake  flow  interactions  are  often  occuring 
simultaneously,  and  the  effects  of  one  may  obscure  the  other.  The 
combined  effects  of  different  plume-disturbing  factors  can  usually  only 
be  determined  by  experiment. 

The  second  wake  flow  examined  in  this  study  is  of  specific 
application  to  the  surface-mining  industry,  but  also  is  of  interest  in 
determining  the  influence  of  terrain  obstacles  on  plumes  in  general. 

On  these  mining  sites,  there  often  exist  large  tailings  ponds  sur¬ 
rounded  by  man-made  dikes,  which  may  become  very  tall  (perhaps  50- 100m) 
over  extended  periods  of  mining.  A  dike  of  this  size  can  easily  be 
the  dominant  topographical  feature  in  the  plant  site  area.  Thus, 
the  wake  flows  created  by  wind  blowing  over  these  dikes  may  significantly 
affect  plume  behaviour.  Investigations  of  dike  effects  are  of  added 
interest,  compared  to  studying  plume  interactions  with  naturally- 
occurring  terrain  obstacles,  because  the  size  and  location  of  the  dikes 
are  subject  to  some  control  by  the  designer. 

Dikes  affect  stack  plumes  mainly  through  the  increased  flow 
turbulence  in  their  wakes.  The  rate  at  which  a  plume  spreads  is 
governed  by  the  amount  of  turbulent  mixing  occurring  in  the  flow  and 
this  is  increased  by  ambient  flow  turbulence.  If  plumes  spread  more 
due  to  dike  wakes,  it  is  likely  that  higher  ground  level  concentrations 
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will  occur,  and  this  effect  would  be  of  importance  to  stack  design. 

Dikes  can  also  be  influential  in  deflecting  ambient  streamlines 
in  the  approach  flow.  In  this  way,  plumes  can  be  shifted  bodily  from 
their  normal  trajectories  by  dike-induced  large  scale  motions  in  the 
atmosphere,  altering  pollution  effects.  It  is  also  possible  that  a 
stack  plume,  trapped  in  the  recirculation  zone  immediately  behind  a 
dike,  could  be  influenced  by  the  increased  mixing  in  this  area. 

Some  simple  analysis  provides  a  hint  of  what  sort  of  effects 
dikes  would  have  on  measured  concentrations.  The  equation 


max 


is  an  estimate  from  the  Gaussian  plume  model  for  the  maximum  value  of 

ground  level  concentration,  C  ,  normalized  by  the  stack  concentration, 

C  .  This  depends  on  the  emission  rate  Q,  the  windspeed  U,  the  plume 

height  h,  and  the  plume  spreads  o  and  a  .  The  maximum  concentration 

y  z 

depends  on  the  ratio  of  o ^  to  a  ,  and  extra  spread  induced  by  dike 

disturbances  should  act  about  the  same  on  a  and  a  ,  causing  their 

y  z 

ratio  to  change  only  slightly.  Thus,  there  should  be  little  effect  on 

the  maximum  gound  level  concentration,  C  .  This  will  later  be 

°max 

shown  from  measurements  of  ground  level  concentrations.  However,  the 

position,  x  ,  at  which  the  maximum  concentration  occurs,  is  very 
v  max 

sensitive  to  plume  disturbances  as  will  also  be  seen. 

The  present  study  employed  physical  modelling  of  these  plume- 
wake  flow  interactions  in  a  wind  tunnel.  The  question,  "Why  a  wind 
tunnel  model?"  can  be  answered  in  part  by  examining  the  disadvantages 
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of  the  major  alternative  methods  for  conducting  such  an  investigation. 

Analytical  methods  for  solving  turbulent  flow  problems  such  as 
this  one  are  not  practical  because  of  the  complexity  of  turbulent 
flow  in  general.  Present  computers  are  orders  of  magnitude  too  small 
to  dynamically  analyse  an  entire  turbulent  flow  including  all  eddy 
sizes  from  the  turbulent  microstructure  to  large  scale  atmosphere  mo¬ 
tions.  Simplifications  can  sometimes  be  used  to  produce  solutions 
that  describe  the  gross  motions  in  a  turbulent  flow,  such  as  assuming 
an  entrainment  velocity  or  eddy  viscosity  model,  but  these  solutions 
are  usually  semi-empirical  in  nature.  Furthermore,  the  geometric 
complexity  of  real  terrain  features,  such  as  buildings  on  a  plant  site 
or  even  a  dike,  tends  to  discourage  attempts  to  obtain  analytical  solu¬ 
tions  . 

Another  alternative  for  solving  the  problems  in  the  present 
study  is  direct  full  scale  measurement.  This  method  has  the  obvious 
advantage  of  realism,  but  also  many  other  large  disadvantages.  Full 
scale  measurements,  as  compared  to  model  measurements,  are  generally 
much  more  expensive  and  difficult  to  perform,  because  of  the  large 
test  areas  involved.  Collecting  full  scale  data  is  also  more  time- 
consuming,  because  physical  processes  such  as  diffusion  are  much  slower 
in  the  atmosphere  than  they  are  in  most  wind  tunnel  models,  and  longer 
sampling  times  are  required.  Also,  in  full  scale  studies,  the  investi¬ 
gator  must  often  wait  for  the  proper  atmospheric  conditions  or  other 
test  configurations  to  be  present.  This  is  not  a  problem  in  wind  tunnel 
studies . 

In  the  atmosphere,  it  is  often  very  difficult  to  test  the  simple 
uniform  flow  cases  that  are  desired  for  comparison  to  analytical  models. 
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For  example,  the  conditions  of  constant  wind  direction  with  time 
and  height  or  uniform  neutrally  stable  atmospheric  stability  are  of 
great  interest  as  baseline  cases  in  the  analysis  of  plume  dispersion. 
However,  these  rarely  occur  in  nature  for  sufficiently  long  periods  of 
time  to  allow  for  these  conditions  to  be  properly  tested.  In  a  wind 
tunnel,  flow  conditions  can  be  closely  controlled  for  as  long  as  nec¬ 
essary,  and  it  is  usually  the  simplest  cases  that  are  the  most  easily 
modelled.  This  latter  feature  can  be  a  disadvantage  also,  because 
sometimes  it  is  desired  to  model  more  complex  flow  conditions  such  as 
positive  atmospheric  stability. 

Wind  tunnel  modelling  of  plumes  has  its  own  unique  problems, 
which  are  mostly  concerned  with  ensuring  an  accurate  simulation  of  full 
scale  processes  despite  certain  modelling  discrepancies.  At  least  as 
much  time  is  usually  spent  in  developing  and  verifying  the  model  as  is 
spent  in  obtaining  the  desired  measurements.  In  this  study,  problems 
were  encountered  with  Reynolds  number  effects  on  the  stack  and  dike 
wake  flows,  because  the  tunnel  was  operated  at  much  lower  Reynolds  num¬ 
ber  than  the  full  scale  atmospheric  flow.  Wake  flows  are  generally 
Reynolds  number  dependent,  and  this  required  a  careful  investigation 
of  the  flow  past  model  stacks  and  dikes,  so  that  their  effect  on  plume 
dispersion  would  be  accurately  simulated. 

The  detailed  objectives  of  the  present  study  were: 

(a)  to  investigate  Reynolds  number  effects  on  flow  past  stacks 
and  dikes  in  a  wind  tunnel  and  ensure  a  correct  model  of 
these  flows 

(b)  to  determine  how  stack  downwash  affects  plume  rise,  and  to 
develop  a  general  rise  model  which  includes  these  effects 
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(c) 


to  determine  how  dikes  affect  plumes  through  measured 
ground  level  concentrations. 


CHAPTER  II 


THE  WIND  TUNNEL  MODEL 

Introduction 

In  this  chapter,  a  description  is  given  of  the  wind  tunnel 
model  and  experimental  technique  used  in  this  study.  An  important 
feature  of  the  model,  its  inability  to  match  mean  flow  Reynolds  num¬ 
ber  to  full  scale,  is  also  discussed. 

Modelling  a  Buoyant  Stack  Plume 

A  hot  gas  emitted  from  a  chimney  has  the  characteristics  of 
both  a  jet  and  a  buoyant  plume,  interacting  with  a  turbulent  cross  flow. 
Modelling  this  very  complicated  flow  requires  first  the  accurate  simu¬ 
lation  of  the  important  turbulent  and  mean  flow  velocity  characteristics 
of  the  atmospheric  boundary  layer  in  which  the  stack  is  immersed. 

Once  this  is  done,  the  approach  flow  to  the  stack  may  be  characterized 
by  one  scaling  velocity,  which  is  usually  the  wind  speed  at  stack  height 

U  . 
s 

Many  different  methods  exist  for  determining  the  correct  model 
values  of  and  the  stack  gas  initial  conditions  to  accurately  model 
a  given  full  scale  stack.  One  of  the  major  differences  among  these 
methods  is  the  way  in  which  the  plume  buoyancy  is  characterized.  Lud¬ 
wig  and  Skinner  (1976)  and  the  present  study  employ  the  flux  modelling 
method,  where  the  flux  of  buoyancy  force  emitted  from  the  stack, 
fg  (P  -  Ps)ws  7rr2J  »  is  considered  to  be  the  relevant  parameter  for 


7 


8 


describing  continuous  buoyant  plumes.  Other  modellers  such  as  Isyumov, 
Jandali,  and  Davenport  (1976)  have  used  the  densiometric  Froude  number, 


U 


Fr  = 


1 

0 

CL 

1 

0 

Q. 

1 _ 

1/2 

g  - ~ -  d 

P 

a 

to  characterize  plume  buoyancy.  The  Froude  number  actually  correctly 
models  buoyancy  effects  for  a  discrete  parcel  of  stack  gas. 

Recent  studies  by  Winkel  (1979)  have  shown  that  the  flux 
modelling  method  accurately  models  plume  trajectories  for  plumes  where 
both  momentum  and  buoyancy  effects  are  important,  over  a  wide  range 
of  plume  density  ratios  p  /p  .  In  the  present  study,  density  ratios 

S  cl 

were  not  varied,  and  the  flux  modelling  method  used  is  summarized 
below. 


The  following  three  parameters  are  equated  between  model  and 
full  scale  (prototype): 

(1)  The  buoyancy  flux  parameter, 


g(pa  ~  ps)d  Ws 

p  U3  x 
a  s 


(2)  The  momentum  flux  parameter, 


2 

ap  w 
s  s 


(3)  The  density  ratio  between  stack  gas  and  the  atmosphere, 

p  /p  .  Simulation  of  full  scale  thermally  induced  den- 
s  a 

sity  ratio  was  obtained  using  molecular  weight-induced 


density  ratio  at  constant  temperature. 
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In  addition,  stack-exit  velocity  profile  shapes  are  modelled 
as  closely  as  possible  to  full  scale  shapes,  characterized  by  their 
kinetic  energy  correction  factors. 


a 


(2-1) 


where  a  is  the  ratio  of  the  mean  kinetic  energy  flux  to  the  kinetic 
energy  flux  found  using  the  profile  mean  velocity,  w  .  In  (2-1), 
wr  is  the  local  stack  gas  velocity  at  a  given  perimeter. 

In  the  present  modelling  method,  all  model  lengths  scale  by 
the  same  factor,  which  is  the  size  reduction  ratio  for  the  model. 

Combining  all  the  above  criteria,  the  modelling  relations  for 
windspeed  Us  and  stack  gas  average  velocity  wg  become: 


U 

U 


SP 

SM 


W 

W 


SP 

SM 


where  the  subscripts 
respectively. 


d„  " 

CM 

r-H 

CL. 

P 

M 

Si 

_v 

d„ 

1/2 

a„ 

P 

M 

d 

an 

M 

P 

1/4 


(2-2) 


3/4 


"M"  and  "P"  refer  to  the  model  and  prototype  values, 


Present  Wind  Tunnel  Model 

The  present  study  was  conducted  in  the  Low  Speed  Wind  Tunnel 
facility  shown  in  Figure  2-1,  located  in  the  Department  of  Mechanical 
Engineering,  University  of  Alberta.  All  experiments  took  place  in  the 
lower  test  section  in  the  configuration  shown  in  Figure  2-2.  The  test 
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section  was  2.44  m  wide  by  1.22  m  high  by  11m  long. 

The  model  plant  site,  shown  in  the  lower  photograph  of  Figure 
2-3,  was  an  800:1  scale  reduction  of  a  typical  oil  sands  strip  mining 
site,  with  flat  terrain  features  except  for  tailings  pond  dikes  tested. 

Two  plume  models,  denoted  as  A  and  B,  were  tested,  with  stack 
conditions  for  each  as  described  in  Table  2-1.  Plume  model  A  used  an 
open  stack  tube  with  laminar  exit  flow  (  a  =  2),  while  model  B  used  a 
stack  tube  with  a  porous  plug  installed  near  the  exit,  which  was 
developed  by  Winkel  (1979).  It  produced  a  turbulent  stack  flow  (a  =1.27), 
with  measured  turbulent  intensities  of  about  i  =  0.2.  Compared  to  full 
scale  stack  flows,  which  are  generally  turbulent  with  flat  velocity 
profiles  (  a  =  1.03),  plume  model  B  was  a  more  accurate  model. 

Reynolds  Number  Mismatch 

In  plume  modelling  with  (2-2),  it  is  not  implied  that  Reynolds 
number,  Re^  =  Ud/v  ,  is  kept  constant  between  model  and  prototype;  in 
fact,  large  mismatches  in  Reynolds  number  often  occur.  These  are 
inherently  present  due  to  the  size  reduction  of  the  model,  but  will 
be  even  greater  if  velocities  in  the  model  are  lower  than  full  scale 
velocities.  For  example,  for  plume  model  B, 


Re 

Re 


dP 

dM 


=  (29.80) (800) (1) 
=  2.38  x  104 


so  that  the  model  Reynolds  number  was  more  than  20,000  times  smaller  than 
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TABLE  2-1:  STACK  CONDITIONS  FOR  MODEL  AND  FULL  SCALE  PLUMES 


Parameter 

Model 

A 

Model  B 

Full  Scale  Equivalent 

Stack  Height, 

22.86 

cm 

22.86  cm 

183m 

Inside  Diameter,  d 

0.99  cm 

0.99  cm 

7.92  m 

Stack  Gas  Density 

Ratio,  P  /p 

5  a. 

0.554 

0.554 

0.554 

Kinetic  Energy 

Factor,  a 

2.0 

1.27 

1.03 

Stack  Gas  Condition 

laminar 

turbulent 

turbulent 

A* 

B* 

Wind  Speeds  at 

0.5 

0.45 

16.7 

13.4 

Stack  Height, 

Us (m/s) 

1.0 

0.90 

33.4 

26.8 

Stack  Exit 

Velocity,  w  (m/s) 

1.006 

0.8604 

46.8 

28.5 

* 


Full  scale  equivalent  from  (2-2)  for  models  A  and  B. 
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the  full  scale  equivalent.  For  plume  model  A,  the  ratio  of  Reynolds 

4 

number  was  2.67  x  10  . 

In  physical  modelling  of  fluid  flows,  a  mismatch  in  Reynolds 
number  means  that  the  model  flow  will  not  be  dynamically  similar  to  the 
prototype.  This  would  certainly  apply  to  the  present  study  considering 
the  large  (fourth  order)  mismatches  calculated  above.  However,  the 
modeller  can  take  steps  to  ensure  that  an  accurate  simulation  is  obtained 
in  spite  of  Reynolds  number  mismatch.  For  the  gross  dynamics  of  the 
plume  motion,  this  was  already  accomplished  by  the  choice  of  the  dis¬ 
torted  model  velocity  scales  using  (2-2) .  For  flows  around  the  plant 
site  obstacles  and  terrain  features  which  could  significantly  affect 
the  intended  measurements  in  the  present  study,  a  careful  investigation 
of  Reynolds  number  effects  was  required,  as  given  in  Chapter  III  for 
stack  cylinders  and  in  Chapter  V  for  dikes. 

Modelling  the  Atmospheric  Boundary  Layer 

Reynolds  number  mismatch  also  had  a  significant  effect  on  the 
method  used  to  simulate  the  atmospheric  boundary  layer  into  which  plumes 
were  emitted. 

The  boundary  layer  was  artificially  created  by  the  use  of  a 
passive  system,  where  the  turbulent  flow  gradually  developed  as  it 
flowed  past  a  series  of  fixed  obstructions.  The  system  used  was  of 
the  type  first  developed  by  Counihan  (1969)  and  consisted  of  the 
following  three  sets  of  turbulence  generating  elements: 


V. 
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(1)  A  row  of  seven  1  m  high  tapered  spikes  extending  to 
near  the  top  of  the  tunnel  section.  These  produced  a 
source  of  vertical-axis  vortex  motions  and  also  a  means 
of  adjusting  lateral  uniformity  of  velocity  across  the 
tunnel  by  varying  the  spike  spacing  and  angle  of  attack. 

(2)  A  low  (10  cm)  barrier  wall  mounted  downstream  of  the 
spikes  and  extending  across  the  width  of  the  tunnel. 

This  provided  a  source  of  horizontal-axis  vortex  motions 
and  a  slowing  of  velocity  near  the  floor  surface  to 
obtain  the  type  of  power  law  mean  velocity  profile  shown 
above. 

(3)  A  10  m  length  of  tunnel  floor  section  covered  with 
roughness  elements,  3  cm  and  6  cm  high.  These  were 
randomly  distributed,  but  in  a  uniform  density  on  each 
roughness  panel,  which  was  120  cm  square.  This  roughness 
provided  further  turbulence  generation  near  the  "ground", 
and  allowed  the  developing  boundary  layer  to  become  later¬ 
ally  and  longitudinally  uniform. 
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Reynolds  number  mismatch  in  the  model  affected  the  design  of 
the  turbulence  generation  system  in  two  ways: 

First,  the  elements  were  made  sharp-edged.  This  meant  that 
flow  separation  points  on  the  elements  were  fixed  at  the  edges,  causing 
their  wake  flows  to  be  very  insensitive  to  Reynolds  number.  Thus,  the 
model  boundary  layer  could  maintain  its  characteristics  over  the 
range  of  model  windspeeds  tested. 

Secondly,  the  elements  were  built  very  large  (24  and  48  m  high 
in  800:1  scale)  in  comparison  to  the  equivalent  full  scale  terrain 
roughness  of  a  forest.  This  exaggerated  flow  disturbance  was  needed 
because  the  boundary  layer  was  to  be  created  with  similar  turbulent 
intensity  as  in  full  scale,  but  at  four  orders  of  magnitude  less  Reynolds 
number.  The  size  comparison  can  be  seen  in  Figure  2-3  where  the  size 
of  the  investigator  is  much  larger  in  comparison  to  the  800:1  plant 
model  shown  in  the  lower  photograph,  than  to  the  spikes,  barrier  and 
roughness  shown  in  the  upper  photograph. 

Turbulence  and  Mean  Velocity  Measurement 

Measurements  of  mean  velocity  and  turbulence  were  taken  with 
linearized  constant  temperature  hot-wire  and  hot-film  anemometers,  as 
shown  in  Figure  2-4.  A  DISA  Type  55D01  hot-wire  anemometer  was  used 
in  the  development  of  the  boundary  layer,  and  a  Thermo-Systems  Inc. 

(TSI)  Model  1050  hot-film  anemometer  was  used  for  all  other  work. 
Single-wire  probes  were  used  to  obtain  mean  windspeed  U  and  turbulent 
intensity  i  =  /  U  ,  while  x-wire  probes  were  used  to  obtain 

velocity  components  in  the  three  coordinate  directions. 

During  the  trial  and  error  adjustments  to  produce  the  simulated 
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boundary  layer,  surface  roughness  elements  were  present  along  the 
entire  10  m  length  of  the  test  section;  later,  the  plant  site  model 
replaced  the  last  5  m  of  roughness  elements,  as  shown  in  Figure  2-2. 

Vertical  profiles  of  mean  and  turbulent  velocities  were  taken 
at  the  stack  location,  and  are  plotted  in  Figure  2-5.  Also,  a  compari¬ 
son  was  made  by  Wilson  (unpublished)  of  characteristic  boundary  layer 
velocity  parameters  between  the  measurements  and  full  scale  data  of 
Counihan  (1975),  as  given  in  Table  2-2.  The  full  scale  condition  in 
Table  2-2  is  representative  of  a  neutrally  stable  atmospheric  boundary 
layer  over  wooded  terrain,  and  the  model  boundary  layer  is  seen  to  be 
a  close  match  to  this  condition. 

Lateral  and  longitudinal  uniformity  of  the  model  boundary  layer 
were  also  tested  and  found  to  be  very  good.  Mean  velocity  was  laterally 
uniform  to  +  6%  of  the  tunnel  centerline  value  over  about  the  central 
80%  of  the  tunnel  width,  while  turbulent  intensity  i  was  laterally 
uniform  to  +_  10%  over  the  central  50%  of  the  width.  Turbulence  intensity 
was  longitudinally  uniform  to  _+  10%  of  the  average  value  along  a  4  m 
length  of  test  section  beginning  near  the  farthest  upwind  point  of  the 
plant  site  location.  This  indicated  that  the  boundary  layer  was  already 
in  equilibrium  before  reaching  the  test  section  where  dispersion  mea¬ 
surements  were  later  carried  out. 

Concentration  Measurement 

The  stack  gases  in  the  model  plumes  were  composed  of  helium- 
air  mixtures.  Helium  served  as  both  the  source  of  buoyancy  and  as  the 
tracer  gas  to  model  contaminant  concentrations.  Conditions  between  the 
model  stack  and  tunnel  were  isothermal  and  isobaric,  therefore  measured 
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TABLE  2-2:  COMPARISON  OF  WIND  TUNNEL  BOUNDARY  LAYER  WITH  FULL 
SCALE  ATMOSPHERIC  DATA  OF  COUNI HAN  (1975) 

(from  Wilson  [unpublished]) 


* 


Data  fron  Hinze  (1975)  p.  643,  p.  729. 
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volume  concentration  ratios  C/Cs  could  be  interpreted  as  mass  concentra¬ 
tion  ratios  X/Xs-  This  is  not  true  for  full  scale  hot  plumes. 

The  procedure  used  to  measure  helium  concentrations,  shown 
schematically  in  Figure  2-6,  was  to  aspirate  a  sample  from  any  desired 
point  within  the  plume  and  feed  it  to  a  heated-element  four  arm  thermal 
conductivity  bridge  of  the  same  type  used  in  commercial  gas-chromatographs. 
The  bridge  was  calibrated  for  helium  volume  concentrations  by  taking 
samples  from  helium-air  mixture  streams  of  measured  volumetric  propor¬ 
tions,  using  rotameters.  The  response  was  linear  at  0.45  millivolts  per 
ppth  of  helium. 

Another  sample,  called  the  "reference",  was  taken  simultaneously 
with  the  plume  sample  from  a  stationary  position  upwind  of  the  plant 
site  model,  to  account  for  background  levels  of  helium  in  the  recir¬ 
culating  tunnel.  The  voltage  output  (imbalance)  of  the  bridge  instan¬ 
taneously  measured  the  net  plume  concentration. 

Using  the  HP5326B  voltage  to  frequency  converter  counting  DVM 
(with  external  clock  timer)  shown  in  Figure  2-6,  four  consecutive  100- 
second  time  averaged  readings  of  net  concentration  were  taken  at  each 
measurement  location.  The  lower  photograph  in  Figure  2-7  shows  the 
concentration  detector  equipment  in  operation,  with  the  traversing 
mechanism  mounted  in  the  tunnel  to  move  the  concentration  sample  probe 
in  a  straight  line  for  profile  measurements. 

An  example  of  a  measured  vertical  concentration  profile  downwind 
of  the  stack  is  given  in  Figure  2-8,  where  the  upper  plot  shows  the 
scatter  of  the  100-second  samples  and  the  lower  plot  shows  the  cumula¬ 
tive  400-second  average  readings. 
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Sampling  Time 

Considering  that  a  plume  travels  downwind  essentially  at  the 
mean  windspeed  U,  the  elapsed  time  t  for  it  to  reach  a  given  downwind 
position  x  is  about  x/U.  The  comparison  of  such  a  time  scale  between 
model  and  full  scale  conditions  gives  an  estimated  comparison  of  how 
fast  physical  processes  are  taking  place  between  the  two  conditions. 
Thus, 


In  the  present  model,  (xp/xM)  =  (dp/dM)  and  Up/UM  may  be  found  from 
(2-2),  giving 


= 

M 

1/4 

y 

tM 

“m 

dM 

1/2 


(2-3) 


Using  (2-3),  the  cumulative  400-second  averaged  readings  of  concentra¬ 
tion,  such  as  shown  in  Figure  2-8,  are  seen  to  correspond  to  estimated 
full  scale  sampling  times  of  2.7  hours  for  plume  model  A,  and  3.0 
hours  for  plume  model  B.  These  estimates  are  probably  too  high,  if 
anything,  because  although  measured  full  scale  plume  spreads  increase 
with  the  amount  of  sampling  time  due  to  wind  direction  variation  with 
time,  measured  wind  tunnel  plume  spreads,  which  are  not  subject  to 
changes  in  weather,  tend  to  become  independent  of  sampling  time  after 
a  certain  period  of  sampling,  and  (2-3)  is  no  longer  applicable. 
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FIGURE  2-1:  LOW  SPEED  WIND  TUNNEL  FACILITY  AT  THE 
UNIVERSITY  OF  ALBERTA 
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FIGURE  2-2:  TEST  SECTION  LAYOUT  AND  DIMENSIONS 
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Turbulence  generating  spikes,  barrier,  and  floor  roughness 


Adjusting  the  Kurz  anemometer  for  windspeed  at  stack  height 


FIGURE  2-3:  TURBULENCE  GENERATION  SYSTEM  AND  PLANT  MODEL 


|0uunj_  puj/w 


22 


si 

o 

1 

CD 


4- 


T 


U) 

c 

cn  ^ 

CD  3 

^  s 

£5 

*6 

< 

J 


FIGURE  2-4:  MEASUREMENT  OF  TURBULENCE  AND  MEAN  VELOCITY 
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Reynolds  Stress,  -u'w',  m2/sec2 


FIGURE  2-5:  TURBULENT  AND  MEAN  VELOCITY  CHARACTERISTICS 

OF  THE  SIMULATED  BOUNDARY  LAYER 
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FIGURE  2-6: 


MEASUREMENT  OF  PLUME  CONCENT AT ION  IN 
THE  WIND  TUNNEL 
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Placing  the  deflector  vane  on  the  dike  crest 


Measuring  a  vertical  profile 


FIGURE  2-7:  DIKE  MODEL  AND  EQUIPMENT  SET-UP 
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CHAPTER  III 


FLOW  PAST  CYLINDRICAL  STACKS 


Introduction 

In  this  chapter,  the  Reynolds  number  sensitivity  of  flow  around 
cylinders  is  examined,  and  applied  to  the  modelling  of  wakes  behind 
cylindrical  stacks.  Because  of  the  mismatch  in  mean  flow  Reynolds 
number  in  the  present  model  (see  Chapter  II),  different  flow  character¬ 
istics  are  expected  for  model  and  full  scale  stack  wakes.  These  dif¬ 
ferences  may  be  reflected  in  a  poor  simulation  of  the  stack  wake  effect 
on  plume  dynamics. 


Cylinder  Pressure  Coefficient 

The  static  pressure  coefficient. 


C 


P  - 


P 

00 


2 

00 


(3-1) 


measures  the  normalized  pressure  at  any  point  in  a  flow  and  is  often 
used  in  pressure  drag  calculations  for  flow  obstacles.  When  applied  to 
the  base  (0  =  180°)  of  a  cylinder,  the  cylinder  base  pressure  coefficient 
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C„,  ,  is  obtained,  where 
Pb 


Pu  -  P 
base  00 


CO 


(3-2) 


Because  the  local  pressure  defect,  (P  -  ) ,  is  fairly  constant  at  all 

angular  positions  on  the  cylinder  past  the  point  of  flow  separation, 

Cpk  effectively  measures  the  magnitude  of  the  pressure  defect  for  the 
entire  initial  wake  region,  and  is  a  comparative  wake  suction  parameter 
for  model  and  full  scale  conditions.  Correct  modelling  of  stack  wake 
effects  on  plume  dynamics  in  a  wind  tunnel  should  then  depend  on  match¬ 
ing  Cpk  to  full  scale,  although  this  is  only  a  simple  criteria  which 
neglects  other  possibly  important  effects  such  as  wake  size  and  turbu¬ 
lence  levels.  Values  of  Cpb  are  either  available  from  previous  studies 
or  are  readily  measurable,  adding  to  the  usefulness  of  this  parameter. 

Infinite  and  Finite  Cylinders 

A  stack  is  a  finite  cylinder,  and  therefore  some  attention  should 
be  paid  to  the  top  end  effect  on  Cpb,  especially  because  the  region  of 
the  stack  wake  flow  near  the  top  of  the  stack  is  in  closest  proximity 
to  the  plume  and  will  exert  the  greatest  influence  on  its  dynamics. 
However,  data  for  infinite  cylinders  was  used  exclusively  in  this  study 
and  assumed  to  apply  to  the  top  portions  of  stacks,  thus  neglecting 
the  end  effect.  This  choice  is  supported  by  the  following  arguments: 

The  end  effect  of  a  stack  top  is  likely  much  less  severe  than 
that  caused  by  a  flat  cylinder  end,  due  to  the  vertical  jet  of  stack 
gas  acting  as  a  cylindrical  flow  obstacle  in  continuation  with  the 
stack.  This  hypothesis  is  strongly  supported  by  the  work  of  Moussa, 
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Trischka,  and  Eskinazi  (1977)  who  conducted  detailed  flow-field  mappings 
of  jets  emitted  from  cylindrical  pipes  in  a  crossflow.  As  shown  in 
Figure  3-1,  they  found  that  "the  leading  surface  of  a  jet,  acting  as 
a  barrier  to  the  crossflow,  is  quantitatively  very  close  to  that  of  a 
rigid  cylinder.  Shed  vortices  are  continuous  with  the  pipe,  and  are 
characteristic  of  the  pipe." 

Although  Moussa's  tests  involved  high-velocity  jets  with  no 
downwash,  all  stack  plumes  would  be  expected  to  exhibit  some  of  the 
above  behaviour  because  even  severely  downwashed  plumes  tested  in  the 
present  study  were  found  to  possess  initial  vertical  rise  of  at  least 
one  diameter  above  the  stack  exit. 

It  therefore  seems  likely  that  flow  near  the  top  of  a  stack 
is  very  similar  to  flow  past  an  infinite  cylinder  and  that  Cp^  values 
are  probably  also  very  close.  Even  if  absolute  values  of  are  not 
exactly  the  same  for  infinite  cylinders  as  for  stack  tops,  a  relative 
comparison  between  values  at  high  full  scale  Reynolds  numbers  and  values 
at  low  model  Reynolds  numbers  in  each  case  would  likely  be  very  similar 
because  the  flows  themselves  are  very  similar.  Thus,  conclusions  re¬ 
garding  Reynolds  number  effects  on  infinite  cylinders  would  be  expected 
to  apply  to  stacks  also. 

Cylinder  Flow  at  Varying  Reynolds  Number 

Variations  in  cylinder  base  pressure  coefficient  with  Reynolds 

number.  Re,  =  Ud  /v  are  closely  tied  to  flow  transitions  that  take 
do 
o 

place  at  various  values  of  Re^  .  It  is  therefore  helpful  to  outline  the 

o 

flow  regimes  that  occur  as  Reynolds  number  increases  from  zero  to  full 
scale  values  for  cylinders. 
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The  first  major  flow  regime  is  the  subcritical  regime,  which 

exists  for  all  Reynolds  numbers  up  to  the  so-called  "lower  transition" 

point,  occurring  in  the  region  2  x  105  <  Re^  <  5  x  105.  In  the  sub- 

o 

critical  regime,  the  flow  separates  from  the  cylinder  while  the  devel¬ 
oping  boundary  layer  is  still  laminar.  This  happens  at  about  0  =  80°, 
causing  a  wake  flow  wider  than  the  cylinder  diameter. 

According  to  the  data  compiled  by  Griffin  (1977),  shown  in 
Figure  3-2  as  dashed  lines,  values  of  base  pressure  coefficient  may  be 

3 

taken  as  approximately  -Cp^  =0.9  below  Re^  =  10  and  -Cpb  =0.8 

3  0 

above  Re^  =  10  .  Data  compiled  by  Sachs  (1972)  shows  that  -C^  increases 
o 

to  a  maximum  value  of  about  ~Cp^  =  1.2  at  the  lower  transition  point. 

In  Figure  3-3,  plots  of  pressure  coefficient,  Cp,  against 

angular  position,  0,  are  given  for  various  flow  regimes,  as  compiled 

by  Roshko  (1960) .  The  subcritical  regime  is  represented  by  the  curve 

from  Fage  and  Falkner  (1931)  at  Re^  =  1.1  x  10^.  The  long  flat  sec- 

o 

tion  beyond  the  separation  point  at  0  =  80°  shows  the  nearly  constant 
value  of  static  pressure  throughout  all  the  rear  surface  area  of  the  cy¬ 
linder  beneath  the  separated  flow.  Values  of  -C^  are  large  in  the 
subcritical  regime  because  the  Cp  curve  flattens  out  due  to  separation 
before  it  has  a  chance  to  return  to  smaller  values,  as  it  did  in  the 
other  curves  in  Figure  3-3. 

The  second  major  flow  regime  for  cylinders  is  the  supercritical 

regime,  and  extends  from  the  lower  transition  point  near  Re^  -  10^  to 

o 

an  upper  transition  point,  which  occurs  within  the  region 

106  <  Re.  <3.5  x  106,  according  to  Roshko  (1960).  In  this  supercritical 
dQ 

regime,  the  cylinder  boundary  layer  undergoes  a  transition  from  laminar 
to  turbulent  flow  prior  to  separation,  allowing  for  a  later  separation 
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point  at  about  0  =  120°.  Also,  -Cp^  drops  steeply  to  a  low  value  of 
about  _cpb  =  0-2,  as  shown  by  Flachsbart's  (1929)  data  in  Figure  3-2. 

The  exact  magnitude  and  location  of  this  drop  are  sensitive  to  approach 
flow  turbulence,  so  that  Flachsbart’s  curve  can  only  be  regarded  as 
typical  of  the  supercritical  regime. 

In  Figure  3-3,  the  data  from  Flachsbart  (1929)  at  Re^  =  6.7  x  10^ 

o 

represents  the  supercritical  regime.  The  small  flat  in  the  curve  near 
0  =  105°  is  probably  evidence  of  a  laminar  separation  "bubble"  [Roshko 
(I960)],  in  which  static  pressure  would  be  expected  to  remain  constant. 
Bursnall  and  Loftin  (1951)  have  shown  that  in  the  supercritical  regime 
for  cylinders,  a  localized  laminar  separation  bubble  occurs  before  the 
turbulent  transition,  reattachment,  and  subsequent  turbulent  separation 
of  the  flow  from  the  cylinder.  Roshko  (1960)  has  observed  that  the 
apparent  disappearance  of  this  bubble  may  be  linked  to  the  upper  transi¬ 
tion  at  about  Re.  =  10^. 

d 

o 

The  third  flow  regime  is  the  transcritical  regime,  and  exists 
at  Reynolds  numbers  above  the  upper  transition  point.  In  this  regime, 
the  turbulent  separation  point  is  observed  by  Roshko  (1960)  to  move 
forward  on  the  cylinder  to  about  0  =  90°.  This  return  to  an  earlier 
separation  and  larger  wake  flow,  as  in  the  subcritical  regime,  is  matched 
by  a  return  in  pressure  coefficient  values  to  subcritical-like  values, 
as  seen  by  Roshko' s  data  in  Figures  3-2  and  3-3.  Roshko  (1960)  also 
suggested  that  no  further  transitions  are  likely  to  take  place  at  higher 
Reynolds  numbers,  because  the  point  of  boundary  layer  transition  to 
turbulence  can  now  only  move  very  gradually  forward  on  the  cylinder. 
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Implications  of  Reynolds  Number  Mismatch 

Because  typical  model  Reynolds  numbers  for  the  present  study 

2  2 

were  in  the  region  3  x  10  <  U^d^/  v  <  8  x  10  ,  while  the  corresponding 
full  scale  Reynolds  numbers  were  in  the  region  8  x  10  <  U  d  /v  <2  x  10, 
it  is  apparent  from  Figure  3-2  that  model  Cpb  values  should  be  very 

close  to  the  full  scale  values  at  about  -CL,  =  0.86. 

Pb 

To  further  test  this  hypothesis,  a  study  was  conducted  in  the 
wind  tunnel  with  a  long  cylinder  immersed  in  the  turbulent  model  boundary 
layer  described  in  Chapter  II.  As  shown  in  Figure  3-4,  the  cylinder 
was  hollow,  with  the  same  outside  diameter,  dQ  =  1.275  cm,  as  the  model 
stack.  However,  unlike  stacks,  it  extended  over  the  entire  test  sec¬ 
tion  height,  to  eliminate  end  effects.  Pressure  taps  were  drilled 
along  the  cylinder  base  at  four  different  heights,  as  shown  in  Figure 
3-4.  Only  one  was  left  open  during  a  specific  test;  the  rest  were  sealed 
with  tape. 

To  obtain  Cpb,  if  was  necessary  to  measure  the  very  small  dif¬ 
ference  between  the  cylinder  base  pressure  and  the  approach  flow  static 
pressure,  which  difference  was  approximately  0.005  cm  of  Ho0  at 
U  =  1  m/sec.  A  Validyne  Model  DP45  differential  pressure  transducer 
was  used,  with  one  input  connected  to  the  inside  of  cylinder  tube, 
which  was  uniformly  at  the  outside  cylinder  base  pressure,  and  the 
other  input  connected  to  a  static  pressure  tube  mounted  in  the  approach 
flow  at  the  same  height  as  the  base  pressure  tap  being  used.  The  trans¬ 
ducer  thus  directly  measured  the  pressure  defect,  p^ase  - 

The  approach  velocity,  U  ,  was  measured  with  either  a  pitot- 
static  tube  (>  2  m/sec)  or  a  Kurz  Model  435  hot  film  anemometer  (0-2 
m/sec),  positioned  at  the  same  height  as  the  pressure  tap  being  used,  so 
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that  the  local  velocity  with  height  was  obtained.  This  local  velocity 
was  clearly  the  correct  velocity  to  use,  as  seen  in  Figure  3-2,  where 
the  measured  values  of  Cp^  vary  smoothly  with  Reynolds  number.  The 
smooth  curve  through  the  data  gives  a  more  accurate  prediction  than  the 
suggested  approximate  values  by  Griffin  (1977).  Turbulence  variation 
with  height  in  the  approach  flow  did  not  significantly  affect  Cp^. 

Because  of  the  fortuitous  match  in  base  pressure-coefficient 
values  between  cylinder  flows  at  model  and  full  scale  Reynolds  numbers, 
model  stack  wakes  are  expected  to  exert  the  same  influence  as  their 
full  scale  counterparts  on  plume  rise,  despite  the  large  Reynolds 
number  mismatch. 
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FIGURE  3-1: 


VORTEX  SHEDDING  FROM  A  JET  [FROM  MOUSSA,  TRISCHKA, 
AND  ESKINAZI  (1977)] 
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FIGURE  3-2:  BASE  PRESSURE  COEFFICIENTS  FOR  INFINITE 
CIRCULAR  CYLINDERS 
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6  (degrees) 


FIGURE  3-3:  PRESSURE  COEFFICIENT  PROFILES  AROUND  CIRCULAR 
CYLINDERS  [FROM  ROSHKO  (I960)] 
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FIGURE  3-4:  MEASUREMENT  OF  CYLINDER  BASE  PRESSURE  COEFFICIENT 


CHAPTER  IV 


PLUME  RISE  AND  STACK  DOWNWASH 


Introduction 

In  this  chapter,  the  stack  wake  effect  on  plume  rise  is  in¬ 
vestigated.  Most  authors,  in  the  study  of  plume  rise,  have  tended 
to  either  ignore  downwash  completely,  or  have  simply  quantified  the 
downwash  effect  on  plume  rise  without  providing  the  means  for  predic¬ 
tion  of  this  phenomenon.  Before  understanding  downwash,  one  must 
first  understand  how  plumes  rise  without  downwash,  and  much  of  this 
chapter  is  devoted  to  this  topic. 

Buoyant  Plume  Rise  in  a  Neutrally  Stable  Atmosphere 

From  full  scale  observations,  Csanady  (1961)  has  divided  the 
plume  trajectory  from  a  stack  into  a  number  of  phases:  an  intial 
vertical  "jet"  phase,  where  momentum  effects  dominate  the  rise,  a 
bent-over  "thermal"  phase,  where  plume  buoyancy  dominates  the  rise, 
a  "break-up"  phase,  where  large  atmospheric  eddies  separate  the 
plume  into  fairly  discrete  puffs,  and  a  "disperse"  phase,  where  smaller - 
scale  atmospheric  turbulence  dominates  the  motions  of  the  dispersed 
parcels  of  stack  gas. 

Although  the  distinction  between  the  first  two  phases  is  use¬ 
ful  in  determining  how  plume  momentum  and  buoyancy  effects  are  dif¬ 
ferent,  it  is  evident  that  the  initial  jet  gradually  and  continuously 
transforms  to  a  bent-over  thermal.  Therefore,  one  continuous  formula¬ 
tion  was  sought  to  describe  plume  rise  from  stack  exit  to  the  termination 
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of  buoyant  plume  behaviour. 

Briggs  (1975)  has  derived  sets  of  governing  differential 
equations  for  the  vertical  and  bent-over  phases  of  rise  for  a  contin¬ 
uous  plume,  from  relationships  for  the  conservation  of  mass,  buoyancy, 
and  momentum.  In  the  case  of  neutral  atmospheric  stability  (no  change 
in  ambient  potential  temperature  with  height) ,  which  is  the  case 
modelled  by  a  uni form -temperature  wind  tunnel,  these  equations  simplify 
greatly,  as  described  below. 

For  both  vertical  and  bent-over  plumes  in  neutral  stability, 

conservation  of  buoyancy  requires  that  the  buoyancy  flux  as  a  function 

of  height,  F  ,  is  constant  and  equal  to  the  source  buoyancy  flux, 
z 

2 

F  =  g([pa  -  P^/Pa  )WSTS  •  The  potential  temperatures  of  both  the  plume 
gas  and  the  ambient  air  are  constant,  due  to  the  assumption  of  adiabatic 
motion  for  the  former,  and  the  condition  of  neutral  stability  for  the 
latter.  Therefore,  the  total  (integrated)  amount  of  buoyancy  force 
exerted  on  the  plume  must  be  constant  with  height.  Adiabatic  plumes 
do  not  include  those  with  condensation  or  evaporation  taking  place,  or 
with  large  radiative  or  chemically- induced  heat  transfers. 

Conservation  of  momentum,  combined  with  the  above  result  from 
conservation  of  buoyancy,  yields  the  same  equation  for  both  vertical 
and  bent-over  plumes: 

d(wV)  =  _L  (4.1} 

dz  — 

w 

In  (4-1),  w  is  the  characteristic  vertical  velocity  of  the  plume  and 
V  is  the  plume  volume  flux  divided  by  tt  (see  Briggs  (1975)  for  exact 
definitions  of  the  integrals  for  the  vertical  and  bent-over  cases) . 
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Equation  (4-1)  was  derived  by  applying  the  Navier-Stokes 
equations  to  a  plume  element  assuming  only  pressure  and  gravity  forces. 
The  viscous  force  term  was  neglected,  bypassing  the  impossible  task 
of  directly  computing  viscous  effects  in  the  turbulent  microstructure 
for  the  entire  flow.  In  order  to  account  for  this  omission,  a  bulk 
assumption  for  the  effect  of  turbulent  mixing  in  entraining  ambient 
air  into  the  plume  was  later  made  in  the  mass  conservation  equations. 
This  is  presently  the  only  practical  way  to  account  for  viscosity  and 
provide  a  closed  set  of  equations.  However,  it  limits  the  analysis 
to  prediction  of  only  the  gross  plume  motion  (centerline  rise). 

The  equations  which  arise  from  conservation  of  mass  take  dif¬ 
ferent  forms  for  vertical  and  bent-over  plumes,  but  have  similar 
character,  as  shown  below. 


dV 

Vertical  Plumes:  dz 

dV 

Bent-over  Plumes:  dz 


1/2 

2  3x(w  V) 

2  32(UV)1/2 


(4-2) 

(4-3) 


In  the  above  equations,  w  and  U  play  similar  roles  in  that 
they  are  the  characteristic  plume  velocity  along  the  plume  axis  in 
each  case.  The  parameters  3-^  and  32  are  the  Taylor  entrainment 
constants,  [Taylor  (1945)]  which  are  related  to  the  bulk  assumption 

for  entrainment,  described  below. 

Consider  a  small  portion  of  a  circular  plume  in  an  arbitrary 

axis  direction,  s,  as  shown  below. 


V 
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The  plume  is  imagined  to  have  a  uniform  velocity  vg  within  a 
characteristic  plume  radius,  r  =  *^v7'T  .  This  is  the  so-called 
’’top-hat"  plume  model. 

Entrainment  of  air  into  the  plume  is  assumed  to  occur  through 
the  sides  at  the  uniform  velocity  v  .  By  the  Boussinesq  approxima- 
tion,  where  the  plume  density  is  taken  as  equal  to  the  atmospheric 
density  in  all  terms  except  those  describing  buoyancy,  mass  conserva¬ 
tion  for  the  plume  element  may  be  evaluated  using  volume  fluxes: 


+  v  2  7T  rds 
e 


A 


+  dv  )ds 
s 


A  +  dA 


With  the  definition  for  volume  "flux"  parameter, 
and  in  the  limit  as  ds  -*•  0, 


V  = 


(1/  tt  ) 


v  ds 
s 


y 


=  2rv 

e 


dV 

ds 


(4-4) 
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Alternately,  (4-4)  can  be  expressed  in  terms  of  plume  radius,  as 


dr 

dt 


(4-5) 


The  Taylor  entrainment  assumption  enters  through  the  choice 

for  v  ,  using  v  =  3,w  '  for  vertical  plumes  and  v  =  $„w  for 

e  e  i  e  2 

bent-over  plumes.  Substituting  these  expressions  into  (4-5),  the 
following  expressions  are  obtained: 

Av 

Non-Buoyant  Vertical  Jets:  —  =  2  ^ 

Buoyant  Bent-Over  Plumes:  -  ^2 

d  z 


The  hypothesis  that  plume  radius  grows  linearly  with  height  has  been 

experimentally  confirmed  for  full  scale  buoyant  plumes  by  Briggs  (1969) 

and  Bringfelt  (1969)  using  photographically  measured  plume  radii. 

Equations  (4-2)  and  (4-3)  are  obtained  by  the  substitution  of 

the  Taylor  assumptions  for  v  into  (4-4) . 

6 

Combined  Rise  Formulation  for  Momentum  and  Buoyant  Rise 

Because  the  conservation  equations  for  momentum  and  buoyancy 
have  been  reduced  to  one  equation,  (4-1),  for  both  vertical  and  bent- 
over  plumes,  it  was  postulated  that  one  continuous  plume  rise  solution 
could  be  obtained  which  would  account  for  both  momentum  and  buoyancy 
effects.  Transforming  (4-1)  into  the  time  domain  with  dt  =  dz/w,  the 
following  is  obtained: 

d(wV)  _ 

~dt"  -  F 


After  integration,  the  solution  is 
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w  V  =  Ft  +  Fm  (4-6) 

The  constant,  F._,  is  the  value  of  the  momentum  flux,  wV,  at  the 
source  (t=0) .  In  this  analysis,  there  is  no  density  included  in  the 
momentum  flux  terms  because  all  plume  densities  are  normalized  by  p  , 

ci 

which  is  the  assumed  constant  density  for  the  flow  by  the  Boussinesq 

approximation.  Thus,  density  in  the  initial  stack  momentum  flux  must 

also  be  normalized,  so  that  Fw  becomes 

M 


The  kinetic  energy  factor,  a  ,  from  (2-1),  is  needed  to  obtain  the 

correct  total  momentum  flux  when  using  the  average  stack  gas  velocity, 

w  . 
s 

A  dilemma  now  appears  in  that  there  are  two  possible  entrainment 

assumptions  in  (4-2)  and  (4-3)  to  combine  with  (4-6) .  Choosing  for  now 

the  model  for  bent-over  rise,  the  plume  radius  follows  dr/dz  =  3^  and  thus, 

r  =  6  Ah  +  r  ,  where  Ah  is  the  plume  rise  above  the  stack  exit  and 
2  o 

r  is  the  initial  plume  radius.  From  the  definition  for  plume  radius, 

2  .  . 

the  volume  flux  is  simply  V  =  Ur  ,  and  the  vertical  plume  velocity  is 
w  =  dz/dt.  Combining  all  of  this,  the  following  is  obtained: 

w  V  =  U  (32  Ah  +  rQ)  (4-8) 

Equating  wV  from  (4-8)  and  (4-6), 

Ft  +  FM  ■  U  '62Ah  +  To)2  ^  (4'9) 
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Assuming  initial  plume  radius  =  0,  (4-9)  may  be  directly  integrated 
to  obtain 


t  + 


(4-10) 


Examining  (4-10), it  is  seen  that  the  first  term  has  precisely 
the  same  form  as  the  "1/3  law"  for  bent-over  jets,  which  has  been 
shown  by  Patrick  (1967)  to  be  accurate  in  predicting  the  rise  of  jets 
in  a  cross-flow.  Also,  the  second  term  has  precisely  the  same  form  as 
the  "2/3  law"  for  buoyant  plumes,  which  has  been  shown  by  Briggs  (1970) 
to  be  applicable  to  the  rise  of  bent-over  buoyant  plumes.  Thus, 

(4-10)  is  very  informative:  it  shows  that  momentum  rise  and  buoyant 
rise  should  add  as  the  sum  of  cubes  in  a  combined  formulation.  It  seems 
logical  then  to  solve  the  dilemma  of  choosing  between  the  two  possible 
entrainment  assumptions  by  modifying  (4-10)  to  include  3^  in  the 
momentum  rise  term,  and  to  retain  3 ^  i-n  buoyancy  rise  term.  With 
this  modification,  and  transforming  the  equation  back  to  x  dependence 
using  x  =  Ut,  the  combined  rise  formulation  is  obtained: 


3F 

J  M  3F  2 

-7-  -7-  x  +  — x 

3^UZ  23p 


(4-11) 


which  is  valid  only  for  a  uniform  velocity  U  with  height.  Equation 
(4-11)  is  no  longer  a  true  solution  of  the  differential  equations,  but 
a  combination  of  the  1/3  and  2/3  laws  for  plume  rise. 

A  finite  initial  plume  radius,  r  ,  can  be  included  in  the  com¬ 
bined  rise  formulation  by  retaining  r^  when  integrating  (4-9) .  The 
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result  is 


Ah  = 


3F 


2  3 

23^ 


3. 


1/3 


r>  (4 


Slawson  (1978)  derived  a  similar  result  including  r^  for  plume  rise  in 

stable  atmospheres.  Equation  (4-12)  was  not  found  to  be  useful  in 

this  study  because  the  correction  term,  became  unreasonably 

large  near  the  stack,  using  r  =  r  .  Therefore,  the  use  of  r  and 

o  s  o 

(4-12)  was  abandoned. 

Equation  (4-11)  may  be  simplified  by  the  definition  of  length 
scales  for  buoyant  and  momentum  rise.  Winkel  (1979)  has  identified 


a  momentum  length,  L^,  where 


Si 


U 


Many  investigators  have  used  a  buoyancy  length,  Lg,  where 


JB 


U 


With  these,  (4-11)  becomes 


Ah' 


3x 

7 


3 


1 


3x 

23: 


B 


Measurement  of  Plume  Rise 

Plume  rise  values.  Ah  =  h  -  h^,  were  measured  for  plume  model 
B  in  the  turbulent  shear  flow  at  windspeeds  =  0.45  and  0.90  m/s. 


-12) 
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Vertical  concentration  profiles,  like  the  one  shown  in  Figure  2-8, 
were  taken  at  a  number  of  downwind  positions  from  the  stack,  from 
x  =  1  cm  (8  m  full  scale)  to  x  =  270  cm  (2.16  km  full  scale),  using  the 
methods  described  in  Chapter  II.  The  measured  profiles  were  fitted  with 
ground-reflected  Gaussian  distributions  (see  Chapter  VI  and  Figure  6-2) 
and  the  plume  height  h  was  defined  as  the  distance  from  ground  level 
to  the  centerline  of  the  unreflected  Gaussian  source  profile  which 
would  generate  the  measured  reflected  profiles,  as  shown  below.  In 
all  cases  in  this  study,  the  plume  height  defined  in  this  way  was 
identical  to  the  height  of  the  point  of  maximum  concentration  on  the 
measured  profiles. 


Measured  plume  heights  for  model  B  are  plotted  in  Figure  4-1. 

Plume  trajectories  were  also  measured  by  Winkel  (1979)  for 
three  variations  on  model  B,  denoted  as  models  Bl,  B2,  and  B3. 

These  plumes  were  tested  with  one  windspeed,  -  0.45  m/s,  and  had 
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the  same  stack  gas  average  velocity  w  ,  density  ratio  p  /p  ,  stack 

S  S3. 

diameter  d,  and  stack  location  as  for  plume  model  B.  The_conditions 
that  were  varied  are  outlined  below: 

(1)  Model  B1  -  The  stack  plug  was  removed,  causing  a  laminar, 
parabolic  stack  velocity  profile  with  a  =  2.0  (a  =  1.27  for 
model  B) . 

(2)  Model  B2  -  The  turbulence  generation  system  and  plant  site 

model  were  removed  from  the  tunnel,  causing  a  uniform  laminar 

crossflow  to  approach  the  stack.  The  stack  height  was  raised 

to  h  =45  cm  (h  =  22.86  cm  for  model  B)  to  avoid  the  shallow 
s  s 

boundary  layer  on  the  tunnel  floor. 

(3)  Model  B3  -  The  turbulence  generation  system  and  plant  site 
model  were  removed  from  the  tunnel  and  replaced  with  a  grid 
of  wooden  slats,  positioned  2.4m  upwind  of  the  stack.  From 
measurements  by  Winkel  (1979)  using  the  TSI  hot-film  anemometer 
system  described  in  Chapter  II,  the  grid  was  found  to  produce 

a  downwind  velocity  field  with  uniform  mean  velocity.  Turbu¬ 
lent  intensity,  i^,  was  vertically  and  laterally  uniform  at 
any  downwind  position  from  the  stack,  but  decayed  longitudinally 

from  about  i  =0.12  at  the  stack  location  to  about  i  =0.05 
u  u 

at  x  =  300  cm.  The  stack  height  was  45  cm. 

A  further  study,  denoted  as  the  laminar  cross-flow  study, 
was  done  using  three  different  laminar  plumes  ( a  =  2)  tested  in  a 
laminar  cross-flow  (bare  tunnel) .  The  inside  stack  diameter  d  was 
1.08  cm  and  the  stack  density  ratio  was  p  /p  =  0.425.  In  this  study, 
the  measured  vertical  concentration  profiles  were  fitted  with  hand- 


:><  . 


48 


sketched  smooth  curves  to  define  the  maximum  concentration  height,  h. 
These  measurements  are  shown  in  Figure  4-2. 

Momentum  Rise  Entrainment  Constant 

Using  data  from  various  investigations,  Briggs  (1975)  has 
shown  that  the  rise  of  momentum  jets  can  be  accurately  predicted  with 
the  1/3  law  when 


(4-13) 


This  semi-empirical  relation  is  physically  reasonable  in  that  jets 
in  a  quiescent  atmosphere  (U  =  0)  should  still  have  a  small  amount 
of  entrainment  due  to  self-generated  mixing,  and  entrainment  should 
increase  with  Us/w  as  higher  windspeeds  cause  increased  disturbance 
to  the  jet  and  force  it  closer  to  the  stack  wake. 

Wilson  (unpublished)  has  suggested  a  more  general  model 
which  is  similar  in  form  to  (4-13),  retaining  the  same  arguments  of 
physical  realism,  but  accounting  for  dependence  of  quiescent  atmos¬ 
phere  jets  on  p  /p  from  the  work  of  Ricou  and  Spalding  (1961). 

s  a 

2  2 

Also,  the  ratio  U  /w  was  replaced  with  <j>M  =  p  w  /p  U  which  was 

S  S  M  S  S  3.  S 

considered  more  relevant  to  plume  rise  in  general.  Wilson's  model 
is 


3 


1 


-1/2 


+  4> 


M 


(4-14) 


In  the  present  study,  which  included  tests  with  variable  a 

and  d>„,  (4-14)  was  found  to  yield  very  accurate  rise  predictions  for 
M 

the  momentum  rise  phase  when  used  with  (4-11).  An  example  of  this  is 
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shown  in  Figure  4-3,  where  measured  plume  rise  data  for  plume  model 
B  at  Us  =  0.45  m/s  is  compared  to  the  prediction  using  (4-11),  shown 
as  a  dashed  line.  In  the  coordinates  used,  the  momentum  rise  portion 
of  (4-11)  is  a  straight  line  with  1/3  slope  and  the  buoyant  rise  por¬ 
tion  of  (4-11)  is  a  2/3  slope  line,  as  shown.  For  this  particular  case, 
momentum  rise  is  seen  to  be  the  dominant  form  of  rise  at  downwind  dis¬ 
tances  less  than  about  x/LD  =  500.  Good  agreement  between  the  data 

D 

and  the  prediction  in  this  region  is  observed,  using  3^  =  0.87  from 
(4-14). 

The  apparent  underprediction  for  the  data  point  near  x/Lg  =  10 
(x  =  1  cm)  is  not  significant,  because  the  absolute  error  of  the  pre¬ 
diction  was  only  about  0.25  cm  in  Ah,  which  was  about  the  resolution 
of  rise  calculations. 

Buoyancy  Rise  Entrainment  Constant  3 2 

Difficulties  were  encountered  in  defining  values  of  3?  to 
be  used  in  the  combined  rise  formulation,  because  of  the  lack  of  a 
simple  comprehensive  model,  equivalent  to  (4-14)  for  jets,  which 
would  predict  32  from  stack  and  atmospheric  conditions. 

The  approach  commonly  taken  in  recent  literature  has  been 
to  assign  a  constant  value  to  32  •  For  example,  in  a  study  using  data 
from  ten  sets  of  buoyant  rise  observations,  Briggs  (1972)  found  an 
overall  best-fit  value  of  39  =  0.6,  using  the  2/3  law.  However,  best- 
fit  values  for  the  individual  studies  ranged  from  =  0.3  to  =  0.9. 
Similar  wide  variations  in  3 2  were  observed  in  the  present  study  when 
matching  measured  rise  with  the  combined  rise  formulation.  Measured 
values  ranged  from  =  0.47  to  =  1.4. 
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The  author  wished  to  avoid  using  (3^  as  a  niere  curve-fitting 
parameter,  because  it  has  physical  significance  in  describing  the  rate 
at  which  a  buoyant  plume  entrains  air  and  grows  with  height,  by 
$2  =  dr/dz.  Therefore,  any  variations  in  $2  should  be  explained  by 
physical  arguments  regarding  entrainment.  For  example,  the  extra 
flow  turbulence  generated  by  a  stack  wake  would  be  expected  to  increase 
entrainment  and  Later,  it  will  be  shown  that  plume  rise  under 

moderate  levels  of  downwash  can  be  predicted  using  adjusted  3 2  values. 

Entrainment  would  also  be  expected  to  be  sensitive  to  wind 
shear  and  turbulence  conditions  in  the  atmosphere,  because  these  are 
related  to  turbulent  mixing.  Also,  initial  stack  flow  conditions  may 
be  important  to  entrainment.  To  investigate  this  effect,  measured 
32  values  from  plume  model  B,  the  laminar  cross-flow  study,  and  models 
Bl,  B2,  and  B3  from  Winkel  (1979)  are  given  in  Table  4-1,  along  with 
stack  and  approach  flow  conditions.  None  of  the  plumes  examined  here 
were  significantly  affected  by  stack  downwash,  so  that  this  effect 
would  not  obscure  the  comparisons. 

From  Table  4-1,  the  effects  of  each  of  the  flow  conditions 
on  32  could  be  seen: 

(1)  Approach  Flow  Shear  -  The  significantly  higher  values  of  32 
for  models  B  and  Bl  than  all  other  values  in  Table  4-1  indi¬ 
cate  that  approach  flow  shear  is  the  largest  single  factor 
in  determining  32*  This  could  be  due  to  two  factors: 
velocity  shear  gradients  enhancing  entrainment  by  moving 
upper  and  lower  plume  layers  past  each  other,  or  higher  mean 
velocities  in  the  upper  regions  depressing  plume  rise  in  a 
way  not  accounted  for  in  (4-11),  which  used  the  windspeed  at 
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stack  height,  U  ,  for  all  calculations. 

(2)  Approach  Flow  Turbulence  -  By  comparing  the  results  for  models 

B2  and  B3  it  is  seen  that  approach  flow  turbulence  alone  does 
not  appear  to  significantly  affect  Entrainment  is  partly 

due  to  buoyancy-induced  self-generated  turbulence,  and  partly 
due  to  atmospheric  turbulence,  and  it  appears  that  the  latter 
is  not  significant  while  the  plume  is  still  rising. 

(3)  Stack  Flow  Conditions  -  By  comparing  the  results  between 
models  B  and  Bl,  or  B2  and  the  laminar  cross-flow  study,  it 
is  evident  that  plumes  with  turbulent  stack  flows  have  higher 

$2  values  than  plumes  with  laminar  stack  flows.  The  differences 
in  a  laminar  cross-flow  are  much  greater  than  those  with  a 
turbulent  approach  flow.  Although  the  stack  profile  shape 
is  accounted  for  in  (4-11)  with  a,  this  does  not  completely 
account  for  the  differences  in  entrainment  when  stack  flows 
vary  widely. 

It  is  clear  that  not  enough  plume  rise  data  was  taken  in  this 
study  to  form  a  comprehensive  model  for  3 2 -  However,  the  results  in 
Table  4-1  allowed  the  author,  in  the  case  of  plumes  which  were  signifi¬ 
cantly  affected  by  stack  downwash,  to  choose  an  appropriate  value  of  32 
and  predict  the  rise  that  likely  would  have  occurred  in  the  absence  of 
stack  wake  effects.  In  this  way,  the  effect  of  stack  downwash  could 
be  quantified. 

Final  Rise 

All  plume  rise  predictions  using  (4-11)  are  only  valid  during 
the  plume  rise  phases  prior  to  plume  break-up  and  domination  by 
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define  a  downwind  position,  xp^na^  at  which  final  rise, 


atmospheric  turbulence.  Thereafter,  plumes  in  turbulent  shear  flows 
are  observed  to  level  off,  or  achieve  final  rise.  It  is  possible  to 

Ah„.  .  occurs, 

Choosing  xp^naj  from  a  measured  trajectory  is  a  somewhat  arbitrary 

decision  due  to  the  continuous  nature  of  rise.  However,  this  did  not 

cause  great  difficulty  in  this  study  because  the  plumes  tested  were 

observed  to  achieve  a  well-defined  final  rise  in  most  cases. 

A  study  was  conducted,  which  involved  the  measuring  of  three 

plume  trajectories  in  the  present  model  boundary  layer.  A  laminar 

stack  flow  condition  was  present,  and  only  one  buoyancy  flux 
4  4  3 

F  =  2.24  x  10  cm  /s  was  used.  The  stack  gas  density  was  constant 
at  p  /p  =  0.784.  By  varying  the  windspeed,  a  fairly  wide  range  of 

S  3. 

final  rise  heights  and  positions  were  tested.  The  results  are  shown 
in  Figure  4-4,  with  smoothed  curves  through  the  measured  rise  values. 

It  was  postulated  that  xfinal  might  depend  on  the  buoyancy 

7 

length,  Ld  =  F/U  ,  since  this  is  an  important  parameter  in  predicting 
B  s 

buoyant  rise.  Choosing  the  simplest  possible  criteria,  xp^na^  “  Lg, 
good  agreement  in  final  rise  position  is  observed  for  the  three 
trajectories  in  Figure  4-4  using 


Xf inal 


2200  Ld 

D 


(4-15) 


This  result  was  in  agreement  with  studies  by  Winkel  (1979)  who  inves¬ 
tigated  final  rise  for  eight  other  plumes,  which  involved  variations  in 

p  /p  as  well  as  a  and  <J>M.  Equation  (4-15)  was  also  found  to  be 
s  a  M 

accurate  for  the  model  B  plume  at  =  0.45  m/sec  as  shown  in  Figure 
4_1.  Thus,  (4-11)  combines  with  (4-15)  to  produce  a  complete  plume 
rise  model  in  neutral  stability  when  downwash  is  not  present. 
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Momentum  Flux  Ratio  (f)^  as  a  Downwash  Parameter 


The  stack-atmospheric  momentum  flux  ratio, 

a  p  w  ^ 
a  s  s 


P  U  ‘ 
a  s 


(4-16) 


has  already  been  used  as  a  plume  modelling  parameter  to  ensure  close 

simulation  of  the  gross  interaction  between  the  stack  jet  and  approach 

flow  in  the  initial  momentum  rise  phase.  Stack  downwash  is  also 

expected  to  be  closely  related  to  (J)^,  because  the  downwashing  pressure 

defect  in  the  stack  wake  is  increased  by  approach  flow  momentum 

1  2 

through  (P,  -  P  )  =  Cn.  p  U  .  Also,  increased  stack  momen- 

^  base  00  2  Pb  a  s 

2 

turn,  ap  w  ,  will  tend  to  lift  the  plume  out  of  the  influence  of  the 
s  s 

stack  wake  and  prevent  downwash.  Because  Cpp)  was  modelled  correctly 
in  the  present  study  (see  Chapter  III),  modelling  <f>^  implies  modelling 
stack  downwash  effects. 

Relating  downwash  to  <f>  is  a  more  general  approach  than  that 
commonly  found  in  most  investigations.  For  example.  Fay,  Escudier, 
and  Hoult  (1969)  used  the  simple  velocity  ratio,  R  =  ws/Us  to  predict 
stack  downwash  effects.  They  found  strong  stack  wake  effects  when 
R  <1.2.  Using  this  as  a  rough  guide,  downwash  might  be  expected  to 
begin  to  be  important  in  about  the  region  1<  <  2. 


Near  Stack  Behaviour  of  Downwashed  Plumes 

Two  sets  of  near-stack  vertical  concentration  profiles  were 

taken  at  varying  values  of  .  These  are  relevant  because  significant 

concentrations  below  the  stack  height  at  small  downwind  distances 

(x  <  2d  )  can  only  be  caused  by  stack  downwash;  plume  spread  is  insigni- 
o 

ficant  at  these  distances. 
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The  first  set  of  profiles  is  shown  in  Figure  4-5.  In  these 
tests,  taken  at  one  diameter  from  the  downwind  stack  edge,  the  constant 
conditions  a  =  2  and  p  /p  =  0.569  were  tested.  Thus,  only  variations 
in  w^  and  were  tested.  In  Figure  4-5,  it  is  seen  that  the  amount 
of  downwash  increased  smoothly  with  decreasing  <j)  .  Downwash  was  very 
severe  at  (j)^  =  0.21  with  measurable  concentrations  five  stack  diameter 
below  stack  height.  Downwash  was  insignificant  at  both  <J>  =  15  and 

<J)^  =3.4  with  no  significant  concentration  measured  below  one  diameter 
below  stack  height.  This  would  indicate  that  when  <j)^  increases  beyond 
a  certain  value,  downwash  virtually  disappears  and  the  value  of  is 
not  important . 

The  second  set  of  profiles  is  shown  in  Figure  4-6.  In  these 
tests,  taken  at  less  than  half  a  diameter  downwind  from  the  trailing 
stack  edge,  a  was  varied  along  with  w^  and  U^.  Again,  decreasing  <f>^ 
increases  stack  downwash,  even  for  plumes  with  different  stack-exit 
velocity  profile  shapes.  Also,  downwash  effects  are  seen  to  be  insigni¬ 
ficant  above  <J>  =  2.57.  It  appears  that  a  small  concentration  will 

always  be  measurable  just  below  the  stack  level.  This  is  probably  due 
to  the  downwashing  of  the  slower-moving  stack  gas  near  the  inner 
wall  of  the  stack. 

Prediction  of  Downwash  Effect  on  Plume  Rise 

Two  plume  trajectories  were  measured  which  displayed  the 

effects  of  stack  downwash. 

In  Figure  4-2,  the  data  at  <f>  =  1.59  in  the  laminar  cross- 

flow  study  shows  a  significant  drop  in  plume  rise  from  the  prediction 
(shown  as  a  dashed  line)  using  (4-11).  Using  &2  =  0.468  and  (3^  from 
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(4-14),  the  other  two  plume  trajectories  in  Figure  4-2  were  accurately 
predicted  by  (4-11),  as  seen  by  the  close  fit  to  the  data.  Except  for 
the  stack  wake  effect,  there  was  no  apparent  reason  why  entrainment 
processes  should  have  differed  for  the  lower  plume.  However,  3^  =  1.0 
was  required  to  fit  the  lower  plume  data,  as  seen  by  the  solid  line 
through  the  data.  This  was  a  moderate  downwash  condition,  where  the 
rise  was  suppressed,  but  could  still  be  modelled  with  (4-11). 

In  Figure  4-1,  the  model  B  plume  at  Us  =  0.9  m/s  (  =  0.643) 

is  clearly  a  case  of  severe  stack  downwash.  The  plume  was  observed  to 
achieve  its  maximum  rise  of  about  Ah  =  1  cm  directly  above  the  stack. 
Thereafter  the  plume  was  entrained  into  the  stack  wake  and  stayed  within 
about  1  cm  of  stack  height. 

The  predicted  rise  of  this  plume,  using  (4-11)  and  (4-15) 
is  shown  for  both  32  =  1-4  and  32  =  0.87.  The  former  was  the  appro¬ 
priate  32  value  for  non-downwashed  plumes  in  the  model  shear  flow,  as 
seen  from  the  upper  plume  at  Us  =  0.45  m/s.  However,  for  low-rising 
plumes,  shear  effects  may  not  be  as  significant,  and  32  =  0.87  is 
the  appropriate  value  from  Table  4-1  for  rise  in  a  non-sheared  turbulent 
flow.  The  difference  between  these  two  predictions  is  seen  to  be 
small  in  this  case. 

Much  more  data  than  that  presented  in  this  study  would  be 

required  to  perform  a  rigorous  correlation  of  downwash  effects  with 

d)  .  However,  from  the  above  results,  downwash  effects  can  be  divided 
M 

into  three  categories: 

(1)  No  Downwash 

For  conditions  of  about  4*^  >  2.0,  stack  wake  effects  are 
minimal,  and  rise  may  be  predicted  from  (4-11)  and  (4-15). 


57 


In  laminar  approach  flows,  final  rise  does  not  occur  and 
(4-15)  is  not  applicable. 

(2)  Rise  Partially  Suppressed  by  Downwash 

For  conditions  of  about  1.0  <(^  <  2.0,  rise  may  be  predicted 

by  (4-11)  using  increased  32  values.  The  following  form  is 

suggested  to  calculate  the  increase  in  3?,  which  is  based  on 

the  assumption  of  no  effect  at  <J>,.  =  2.0  and  the  data  in 

M 

Figure  4-2. 

32  (with  downwash)  =  32  (6  -  2.5  <|>M) 

for  1.0  <  (ft,  <  2.0 
M 

This  is  only  a  rough  estimate,  but  has  the  correct  trend 
with 

M 

(3)  Rise  Completely  Suppressed  by  Downwash 

For  conditions  of  about  <J>,,  <  1.0,  stack  wake  effects  will  be 

M 

severe,  and  a  reasonable  rise  estimate  is  Ah  =  0. 

With  the  above  downwash  categories,  and  the  rise  models  (4-11) 
and  (4-15),  rise  may  be  accurately  predicted  for  plumes  with  both 
momentum  and  buoyancy  in  neutrally  stable  atmospheres. 


Buoyancy  Flux  F  =  0.92185  x  104  cm4/sec3 
Kinetic  Energy  Factor  a  =  1.27 
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FIGURE  4-1:  MEASURED  AND  PREDICTED  PLUME  TRAJECTORIES  FOR 
MODEL  B  PLUMES 
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FIGURE  4-2:  MEASURED  AND  PREDICTED  PLUME  TRAJECTORIES  FROM 
THE  LAMINAR  CROSS-FLOW  STUDY 
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COMBINED  RISE  FORMULATION 
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FIGURE  4-4:  FINAL  RISE  OBSERVATIONS  IN  THE  WIND  TUNNEL 
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FIGURE  4-5:  NEAR-STACK  VERTICAL  CONCENTRATION  PROFILES  AT 
x  =  1  5 

'  dQ  FOR  PLUMES  UNDER  THE  INFLUENCE  OF 
STACK  DOWN WASH 
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CHAPTER  V 


FLOW  RECIRCULATION  BEHIND  DIKES 


Introduction 

It  was  seen  in  Chapter.  Ill  that  the  Reynolds  number  is  signifi¬ 
cant  in  determining  the' 'proper  modelling  of  stack  wake  suction  and  plume 
downwash.  In  this  chapter  we  focus  our  attention  on  the  wake  flows 
behind  much  larger  two-dimensional  flow  obstacles:  the  dikes  which 
surround  tailings  ponds  typically  occurring  on  oil  sands  strip-mining 
sites . 

Dikes  are  similar  to  stack  cylinders  in  that  the  large  mismatch 
in  Reynolds  number  between  the  wind  tunnel  model  and  full  scale  is  again 
found  to  be  a  significant  potential  source  of  error.  However,  unlike 
cylinder  flow,  the  flow  around  such  arbitrarily-shaped  obstacles  as 
dikes  is  more  difficult  to  predict.  It  was  found  that  these 
flows  require  modification  to  account  for  Reynolds  number  effects  before 
the  dikes  can  be  used  in  a  wind  tunnel  dispersion  model  without  serious 
modelling  error. 

Previous  Investigations  of  Recirculation 

The  results  of  several  previous  studies  of  recirculating  flows 
are  helpful  in  demonstrating  the  following: 

(a)  The  amount  of  sensitivity  of  recirculating  flows  to  Reynolds 
number  for  different  obstacle  shapes. 

(b)  Problems  with  quantifying  the  size  and  influence  of  a  recircu¬ 
lation  zone. 
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(c)  Factors  other  than  Reynolds  number  that  may  require  modelling 
for  proper  simulation  of  a  dike  wake. 

For  sharp-edged  objects,  the  wake  flow  is  quite  insensitive  to 

Reynolds  number  because  separation  always  occurs  at  the  leading  sharp 

edge,  thus  essentially  fixing  the  size  of  the  resulting  recirculation 

zone.  The  collected  data  in  Table  5-1  for  measured  recirculation  lengths 

L  /H  behind  thin  sharp-edged  walls  supports  this  claim.  L  is  the  dis- 
c  c 

tance  from  the  wall  to  the  point  of  reattachment  of  the  recirculating 
zone  on  the  downstream  surface,  and  H  is  the  height  of  the  wall.  It  is 
seen  in  Table  5-1  that  the  decreasing  trend  of  L^/H  down  the  table  does 
not  appear  to  be  correlated  with  changes  in  Re^. 

At  the  other  extreme,  gently  curved  obstacles  will  also  be 
Reynolds  number  insensitive  because  separation  and  recirculation  will 
not  normally  occur.  However,  many  commonly-occurring  flow  obstacles 
which  are  rounded  but  not  sharp-edged  exhibit  dramatic  sensitivity  to 
Reynolds  number,  as  discussed  in  Chapter  III  for  circular  cylinders.  The 
dikes  tested  in  this  study  also  fall  into  this  latter  category,  requiring 
careful  investigation  of  the  effect  of  Reynolds  number. 

In  attempting  to  measure  the  reattachment  length  Lc  behind  an 
obstacle,  problems  with  flow  fluctuation  are  often  encountered.  Instan¬ 
taneously  observed  reattachment  lengths  may  vary  considerably  around  a 
mean  value,  as  found  by  Wilson,  Winkel,  and  Neiman  (1979)  in  a  flow 
visualization  study  of  recirculation  zones  behind  model  walls  and  dikes 
in  a  water  channel.  Figure  5-1  shows  the  results  of  two  such  tests  where 
the  fluctuating  reattachment  zone  was  found  by  injecting  small  puffs  of 
dye  near  the  floor  surface  and  observing  whether  the  puffs  were  carried 
upstream  or  downstream.  The  two  positions  beyond  which  the  puffs  always 
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TABLE  5-1:  FLOW  REATTACHMENT  LENGTH  BEHIND  TWO-DIMENSIONAL  WALLS 
IMMERSED  IN  A  BOUNDARY  LAYER 


Investigator 

R 

Re..  =  - 

H  v 

6 

H 

L 

c 

H~ 

Arie  and  Rouse 
(1956) 

4 

3  x  10 

0 

16.8 

Good  and  Joubert 
(1968) 

1.76  x  105 

0.75 

13.3 

Chang  (1966) 

4.7  x  104 

6 

12.8 

Plate  and  Lin 
(1964) 

4 

2  x  10 

6 

12 

Wilson,  Winkel 
and  Neiman  (1979) 

7.3  x  103 

5 

9.5 

Present  study 

6.3  x  104 

8.6 

9.05 
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moved  upstream  or  downstream  defined  the  boundaries  of  the  zone.  From 
Figure  5-1  this  zone  was  typically  30-40%  the  length  of  the_mean  cavity 
length  L^.  The  implication  of  this  is  that  simple  flow  visualization 
techniques  for  measuring  by  observing  how  far  downstream  recirculation 
occurs  may  overpredict  by  perhaps  as  much  as  20%.  A  far  downstream 
fluctuation  of  will  always  leave  dye  behind  when  it  recedes,  and  this 
may  obscure  the  true  mean  position.  These  results  led  the  author  to 
exercise  caution  when  comparing  measured  reattachment  lengths  from  dif¬ 
ferent  experimental  techniques,  and  pointed  to  the  need  for  a  technique 
to  measure  the  true  mean  length. 

Recirculating  flows  may  also  be  sensitive  to  approach  flow 
turbulence;  greater  turbulence  levels  being  expected  to  induce  earlier 
reattachment.  The  parameter  5/H  for  a  recirculation  study,  where  6  is 
the  boundary  layer  thickness,  can  serve  as  a  measure  of  this  effect  since 
turbulence  intensity  decays  with  height  in  a  boundary  layer  to  small 
values  at  height  6.  Thus,  if  6/H  is  large,  the  recirculating  flow  will 
experience  high  turbulence  levels  and  should  have  a  smaller  value  of 
L  /H.  This  general  trend  is  seen  to  occur  for  the  data  shown  in  Table 
5-1,  indicating  the  need  to  properly  model  the  parameter  6/H  and  the 
corresponding  boundary  layer  turbulence  levels  for  the  recirculation 
tests  on  model  dikes  in  the  present  study. 

Wind  Tunnel  Recirculation  Tests 

Tests  were  conducted  in  the  wind  tunnel  boundary  layer  and  plant 
site  model  described  in  Chapter  II,  Three  dike  sizes,  referred  to  through¬ 
out  this  study  as  large,  medium,  and  small,  were  tested,  with  dimensions 
given  in  Table  5-2.  The  dikes  were  always  mounted  with  their  crests 


TABLE  5-2:  MODEL  DIKES  TESTED 
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perpendicular  to  the  flow  in  the  "upstream"  dike  position  shown  in 
Figure  2-  2  .  A  10.4  cm  high  sharp-edged  wall  was  also  tested  for  recir¬ 
culation. 

The  dikes  were  always  tested  with  the  same  base  surface  level 
on  both  sides,  which  implies  a  low  water  level  in  the  tailings  pond. 

Flow  Reattachment  Detector 

Because  of  fluctuation  of  the  reattachment  point  shown  in  Figure 
5-1,  a  technique  was  developed  to  measure  the  true  mean  reattachment 
length  of  a  recirculating  flow.  This  consisted  of  a  probe  mounted  on  the 
floor  surface  downwind  of  the  obstacle,  and  an  external  measurement  and 
monitoring  system,  shown  schematically  in  Figure  5-2.  Its  principle  of 
operation  was  to  inject  helium  tracer  gas  from  a  central  source  tube 
perpendicularly  onto  the  surface  and  to  allow  the  local  surface  flow  to 
carry  the  tracer  to  either  one  or  both  of  two  helium  detector  probes 
positioned  1  cm  upwind  and  downwind  of  the  source.  The  two  detector 
probes  were  simply  tubes  into  which  samples  were  simultaneously  aspir¬ 
ated  and  fed  into  the  sample  and  reference  sides  of  the  helium  detector 
bridge  described  in  Chapter  II.  The  mean  direction  of  surface  air  flow 
(upwind  or  downwind)  was  indicated  by  the  probe  which  gave  the  largest 
reading  over  a  suitably  long  averaging  time.  Actually,  the  two  detec¬ 
tor  readings  were  instantaneously  subtracted  so  that  the  polarity  of 
the  net  reading  yielded  the  flow  direction.  When  the  detector  was  lo¬ 
cated  so  that  the  net  reading  was  zero,  this  position  was  taken  as  the 
mean  point  of  reattachment  for  the  recirculation  zone.  Example  readings 
for  two  flow  conditions  are  given  in  Figure  5-2  for  measurement  locations 
A  and  B.  In  the  first  example  with  ReH  =  1990,  the  mean  point  of 
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reattachment  must  lie  between  A  and  B  because  the  surface  air  flow  was 
upwind  at  A  and  downwind  at  B,  as  shown  in  the  figure.  In  the  second 
example  with  Re^  =  14200,  the  recirculation  zone  was  much  smaller,  with 
downwind  flow  at  both  A  and  B. 

By  moving  the  probe  and  observing  the  readings  as  above,  the 
reattachment  length  could  be  measured  for  a  given  dike  and  approach 
velocity.  For  greater  sensitivity  in  reading,  the  helium  detector  output 
was  usually  given  a  DC  offset  and  measured  on  a  voltage  to  frequency 
averaging  digital  voltmeter  for  100  sec.  samples. 

The  detector  system  had  one  operating  disadvantage  in  that  it 
could  normally  only  be  used  on  a  horizontal  surface.  On  an  inclined 
surface  such  as  a  dike  face,  the  buoyant  helium  tracer  gas  tended  to 
migrate  uphill  and  biased  the  reading  significantly.  This  was  not  a 
serious  problem  because  at  most  Reynolds  numbers  the  mean  reattachment 
point  extended  beyond  the  downwind  edge  of  the  dike  base. 

Experimental  Results 

Using  the  flow  reattachment  detector,  reattachment  lengths  were 
measured  for  the  three  dikes  over  a  range  of  Reynolds  numbers  from 
lower  than  typical  model  values  to  as  high  as  the  tunnel  would  safely 
permit.  The  results  are  plotted  in  Figure  5-3  for  two  different  choices 
for  Reynolds  number,  along  with  two  measurements  in  a  water  channel  from 
Wilson,  Winkel,  and  Neiman  (1979). 

uhh 

The  first  choice  for  a  correlation  parameter,  Re^  =  — ,  was 
made  using  the  dike  and  tunnel  parameters  that  were  varied  during  the 
tests.  It  is  seen  from  Figure  5-3  that  this  yields  only  a  fair  correla¬ 
tion  of  the  data.  This  seemed  to  indicate  that  the  proper  Reynolds 
number  to  correlate  reattachment  lengths  should  not  depend  on  the  obstacle 
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itself,  but  only  on  the  approach  flow.  The  necessary  information  on 

the  dike  height  would  then  be  contained  only  in  L  /H.  The  second  choice. 

U^Q  c 

ReQ  =  used  the  momentum  thickness  of  the  approach  flow  boundary 

layer  0  =  13.14  cm,  found  from  integration  of  the  wind  tunnel  vertical 
velocity  profile  given  in  Chapter  2.  The  boundary  layer  thickness 
6  =  94  cm  was  used  to  calculate  from  the  power-law  fits  to  the  velo¬ 
city  profile,  starting  with  the  measured  values  of  Uu  in  each  case.  In 
Figure  5-3  it  is  seen  that  ReQ  correlates  the  data  much  better  than  RetI, 
although  it  is  impossible  to  say  whether  some  other  approach  flow  para¬ 
meters  than  0  and  would  provide  the  best  correlation  in  general, 
because  only  one  boundary  layer  was  tested. 

Two  explanations  aTe offered  for  the  20-30%  increase  in  L  /H  for 

c 

the  water  channel  data  and  the  large  dike  data  over  the  curve  in  Figure 
5-3  for  the  two  smaller  dikes. 

One  possible  explanation  involves  the  effect  of  approach  flow 
turbulence,  which  tends  to  retard  reattachment.  This  is  analogous  to 
the  arguments  used  to  explain  the  variations  in  the  data  of  Table  5-1. 
Clearly,  6/H  increases  as  L^/H  decreases  for  the  data  in  Figure  5-3  and 
this  may  be  due  to  a  decrease  in  turbulence  as  higher  dikes  encounter 
lower  levels  of  boundary  layer  turbulence. 

The  second  explanation  involves  the  effect  of  tunnel  area  block¬ 
age  by  the  dikes.  Any  obstacle  in  a  fixed-area  tunnel  reduces  the  effec¬ 
tive  flow  area  and  causes  an  acceleration  past  the  obstacle  with  a 
corresponding  static  pressure  drop.  When  the  flow  decelerates  past  the 
obstacle  the  pressure  returns  to  the  higher  ambient  levels.  These  condi¬ 
tions  of  deceleration  against  an  adverse  pressure  gradient  are  exactly 
those  that  enhance  separation  and  retard  reattachment.  This  effect  was 
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measured  by  Good  and  Joubert  (1968)  who  found  that  rear  separation  bubbles 
behind  walls  in  turbulent  boundary  layers  were  affected  by  ^pressure 
gradients  induced  by  tunnel  blockage.  Interestingly,  for  fixed  area 
tunnels  where  6  is  close  to  the  tunnel  height,  such  as  in  the  present 
study  and  the  water  channel  study,  5/H  and  tunnel  blockage  are  inversely 
proportional  to  each  other.  Therefore,  both  turbulence  levels  and  block¬ 
age  varied  with  5/H,  so  it  was  not  possible  to  assign  the  effects  on 
reattachment  to  one  or  the  other,  although  tunnel  blockage  is  probably 
the  more  likely  explanation. 

Figure  5-4  shows  the  correlation  of  L^/H  with  tunnel  area  block¬ 
age  as  taken  from  the  data  in  Figure  5-3  at  ReQ  =  9600.  From  this  it  may 
be  concluded  that  recirculation  will  not  be  affected  by  tunnel  blockages 
up  to  about  5%,  but  may  increase  significantly  when  the  blockage  is 
greater. 

Implications  of  Reynolds  Number  Mismatch 

Figure  5-3  shows  that  the  wake  flows  behind  the  dikes  in  the 

present  study  are  very  sensitive  to  Reynolds  number,  as  well  as  second- 

order  effects  such  as  tunnel  blockage.  Recirculation  zones  appear  to 

shrink  to  negligible  size  when  Reynolds  numbers  increase  to  about 

Re  >  2  x  105,  one  to  two  orders  of  magnitude  above  model  values,  which 
9 

were  typically  from  ReQ  =  10^  to  10.  Because  the  corresponding  full 

7 

scale  values  are  of  about  order  10  ,  it  follows  that  the  present  dike 
model,  if  used,  would  greatly  exaggerate  the  size  and  influence  of  the 
dike  recirculation  zone.  Full  scale  measurements  by  Eliseev  (19 / o) 
support  this,  indicating  that  for  full  scale  hills  about  10  steeper 
than  the  present  model  dikes,  L^/H  =  2  for  Re^  =  10  .  This  means  that 
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full  scale  dikes  would  be  expected  to  have  Lc/H  less  than  2  at  the  same 
Reynolds  number.  _ 

Because  it  was  the  purpose  of  this  study  to  gauge  the  effect  of 
dike  recirculation  zones  on  air  pollution  dispersion,  this  modelling 
error  was  considered  very  serious  and  could  not  be  tolerated.  Model 
Reynolds  number  could  not  be  altered  due  to  stack  plume  modelling  con¬ 
siderations,  so  a  remedy  was  sought  by  modifying  the  dike  flows  as  they 
existed. 

Dike  Deflector  Vane 

A  thin  metal  deflector  vane  was  mounted  along  the  entire 
length  of  the  model  dike  crests  to  provide  flow  in  the  model  close  to 
that  of  full  scale  conditions.  Since  elimination  of  the  large  recircu¬ 
lation  zones  was  the  goal,  the  vane  was  designed  to  force  the  flow  down 
the  leeward  side  of  the  dikes  as  shown  below. 


DIKE  AND  VANE  CROSS-SECTION 


74 


Flow  visualization  studies  were  undertaken  in  the  wind  tunnel 
to  qualitatively  evaluate  and  optimize  the  performance  of  the  vane.  These 
were  done  using  titanium  tetrachloride  smoke  released  upstream  of  the 
dikes.  Only  the  large  and  medium  dikes  were  tested  (or  used  again  in 
this  study),  and  the  same  vane  was  used  for  both  dikes.  After  finding 
the  optimum  spacing  between  the  dike  crest  and  vane,  the  vane  was  ob¬ 
served  to  almost  completely  eliminate  recirculation  zones  behind  both 
dikes  at  all  model  Reynolds  numbers.  Some  typical  results  are  given  in 
Figure  5-5  where  recirculation  zones  were  sketched  from  flow  visualiza¬ 
tion.  At  an  intermediate  Reynolds  number  between  those  in  the  top  two 
sketches  where  flow  recirculation  zones  were  both  large,  the  bottom 
sketch  with  the  vane  in  use  shows  recirculation  isolated  to  two  small 
areas:  one  near  the  dike  base,  and  one  immediately  behind  the  vane. 

These  zones  diminished  with  increasing  Reynolds  number  and  were  never 
observed  to  be  larger  than  those  shown  in  Figure  5-5. 

Use  of  the  helium  flow  reattachment  detector  confirmed  the  suc¬ 
cess  of  the  vane.  Tests  showed  that  net  upwind  surface  flow  could  not 
be  measured  at  a  distance  5  cm  downwind  of  the  large  dike  base  at  any 
tunnel  speed.  Also,  at  distances  of  2.5  cm  behind  the  medium  dike  base 

and  1.3  cm  behind  the  large  dike  base,  recirculation  was  measurable 

3 

only  when  Reynolds  number  was  reduced  to  Reg  <  4  x  10  . 

Comparison  of  Dike  Model  to  Full  Scale 

All  of  the  above  tests  for  dike  models  considered  only  the 
size  of  the  recirculation  zone  as  an  initial  gauge  of  the  validity  of 
the  model.  To  further  verify  that  the  vane  was  operating  correctly, 
the  effects  of  the  model  dikes  on  downstream  mean  and  turbulent  flow 


* 


75 

patterns  were  measured  and  compared  to  the  expected  full  scale  effects. 
In  the  absence  of  full  scale  data,  the  substitute  used  was  to  operate 
the  tunnel  at  a  high  Reynolds  number,  Re^  =  1.3  x  10^,  with  the  vane 
removed.  Figure  5-3  shows  that  the  recirculation  zone  behind  the 
dikes  has  already  decreased  to  negligible  size  at  this  Reynolds  number, 
indicating  a  proper  simulation  of  full  scale  conditions.  This  data 
will  be  referred  to  below  as  the  "full  scale  condition".  The  typical 

model  condition  tested  for  comparison  was  operated  with  Re„  =  1.1  x  10^ 

y 

(U„  =  0.72  m/s)  using  the  dike  with  and  without  the  vane. 

With  a  linearized  single  hot-film  anemometer,  vertical  profiles 
of  mean  velocity  U  and  turbulent  intensity  i  =  17  /  U  were  taken 
at  a  location  x  =  7H(82.25  cm)  downstream  of  the  large  dike,  as  shown 
in  Figure  2-  2.  The  equipment  set-up  for  turbulence  measurement  is 
shown  in  Figure  2-  4. 

In  Figure  5-6  the  results  are  presented  for  mean  velocity  pro¬ 
files  at  the  "full  scale"  condition,  with  and  without  the  dike  present. 
The  normalizing  mean  velocity  Ure£  was  taken  as  the  velocity  at  Z  =  5H 
for  each  case,  well  above  the  immediate  flow  effects  of  the  dike.  It 
is  seen  from  Figure  5-6  that  at  "full  scale"  conditions,  the  dike  pro¬ 
duces  very  little  disturbance  of  the  mean  velocity  profile  at  x  =  7H. 

A  good  model  would  therefore  be  expected  to  exhibit  similar  behaviour. 

Figure  5-7  shows  the  normalized  mean  velocity  profiles  for 
the  two  model  dikes  with  and  without  the  vane,  compared  to  model-speed 
no-dike  case.  It  is  apparent  that  the  model  dike  with  the  vane  produces 
a  close  match  in  velocity  profile  to  the  no-dike  case  while  the  model 
dike  without  the  vane  yields  a  very  poor  match.  This  can  also  be  seen 
by  examining  power  fits  to  the  above  profiles  as  given  in  Table  5-3. 
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TABLE  5-3:  VELOCITY  PROFILE  POWER  LAW  FITS  BEHIND  DIKES  AT  X  =  7H 


U 

U 


ref 


n 


Condition 

Dike 

Height  Range  (cm) 

n 

"Full  Scale" 

Re0  =  1.3  x  105  ■ 

Large  (no  vane) 

0-10 

10-50 

0.094 

0.240 

No  Dike 

0-10 

10-50 

0.102 

0.230 

Model 

Re0  =  1.1  x  10^ 

Large  (no  vane) 

0-20 

20-40 

0.33 

Large  (with  vane) 

0-15 

15-40 

0.205 

0.300 

No  Dike 

0-20 

20-80 

0.182 

0.309 
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The  model  dike  with  vane  and  the  model  no-dike  case  match  closely  in 
power  law  as  is  the  case  in  "full  scale",  while  the  model  dike  case 
without  the  vane  does  not  yield  to  a  power  law  fit  at  all,  being  strongly 
linear  in  velocity  with  height  up  to  about  two  dike  heights. 

In  Figure  5-8  the  vertical  turbulence  profiles  for  the  "full 
scale"  and  two  model  cases  are  plotted  for  comparison.  The  turbulent 
disturbance  parameter  used  is: 

Ai  i  (with  dike)  -  i  (without  dike) 

u  _  _u _  u  _ _ 

1U  i^(without  dike) 

This  parameter  is  the  fractional  increase  in  intensity  due  to  the  dike 
wake,  or  in  other  words,  the  amount  of  turbulence  added  to  the  flow  by 
the  dike.  One  would  expect  this  ratio  to  be  closely  matched  to  full 
scale  values  by  a  good  model.  Figure  5-8  shows  that  the  model  dike 
with  vane  does  produce  a  close  match  to  the  "full  scale"  case  while  the 
model  dike  without  the  vane  adds  far  too  much  turbulence  to  the  lower 
levels  of  the  flow,  as  would  be  expected  from  the  measured  recirculation 
zones  produced  by  this  model. 

Summary 

The  mismatch  in  Reynolds  number  between  model  and  full  scale 
gives  model  dikes  with  grossly  exaggerated  recirculation  zones.  The 
use  of  a  deflector  vane  mounted  on  the  model  dike  crests  corrects  this 
problem,  producing  a  dike  flow  which  is  a  close  match  to  expected  full 
scale  in  downstream  mean  velocity  and  turbulence  profiles.  The  model 
dike  with  vane  can  be  used  to  accurately  simulate  the  effects  of  dikes 


on  atmospheric  disperson. 
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FIGURE  5-1:  FLOW  SEPARATION  REGIONS  OBSERVED  IN  A  WATER 

CHANNEL  [FROM  WILSON,  WINKEL  AND  NEIMAN  (1979)] 
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FIGURE  5-2:  MEASUREMENT  OF  REATTACHMENT  LENGTH  WITH  THE  HELIUM 
FLOW  DIRECTION  DETECTOR  [FROM  WILSON,  WINKEL  AND 
NEIMAN  (1979)] 
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Hill  Height  Reynolds  Number  in  Approach  Flow 
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FIGURE  5-3: 


REYNOLDS  NUMBER  EFFECTS  ON  FLOW  REATTACHMENT  BEHIND 
MODEL  DIKES  [FROM  WILSON,  WINKEL  AND  NEIMAN  (1979)] 
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FIGURE  5-4:  TUNNEL  BLOCKAGE  EFFECTS  ON  REATTACHMENT  LENGTHS 
BEHIND  MODEL  DIKES 
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Ree  =  2.7  x  103 


Ree  =  6.7  x  103 


Re^  =  5.4  x  103 
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FIGURE  5-5:  EFFECT  OF  THE  DIKE  DEFLECTOR  VANE  ON  DIKE  WAKE 
FLOW  FROM  FLOW  VISUALIZATION  IN  THE  WIND  TUNNEL 
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SURE  5-6: 


VERTICAL  MEAN  VELOCITY  PROFILES  BEHIND  A  MODEL  DIKE 
AT  "FULL  SCALE"  CONDITIONS 


Model  Height  Above  Ground  (cm) 
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FIGURE  5-7:  VERTICAL  MEAN  VELOCITY  PROFILES  BEHIND  MODEL  DIKES 
AT  TYPICAL  MODEL  CONDITIONS 
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Model  Height  Above  Ground  (cm) 
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FIGURE  5-8:  VERTICAL  PROFILES  OF  DIKE  DISTURBANCE  ON  TURBULENT 
INTENSITY  BEHIND  MODEL  DIKES 


CHAPTER  VI 


INTERACTION  OF  DIKES  WITH  STACK  PLUMES 


Introduction 

In  this  chapter  the  interaction  between  dike  wakes  and  the 
disperson  of  stack  plumes  will  be  investigated  by  examining  the  effects 
of  dikes  on  ground-level  concentrations  (GLC) .  Increases  in  GLC 
due  to  any  factor  are  of  prime  interest  to  the  stack  designer  because 
of  the  emphasis  on  maximum-allowable  GLC  levels  in  air  pollution  legis¬ 
lation.  A  simple  semi-empirical  theory  will  be  developed  for  estima¬ 
ting  dike  effects  on  GLC  values. 

Experimental  Procedure 

Profiles  of  GLC  were  taken  along  a  straight  line  downwind 
of  the  model  stack  shown  in  Figure  2-2.  Various  dike  sizes  and  con¬ 
figurations,  including  dikes  placed  upstream  and  downstream  of  the 
stack  were  tested,  as  summarized  in  Table  6-1.  The  dikes  were  always 
oriented  perpendicular  to  the  flow  and  tested  with  the  dike  flow  de¬ 
flector  vane  described  in  Chapter  V.  Plume  models  A  and  B  from  Table 
2-1  were  used,  where  model  A  had  a  laminar  stack  gas  and  model  B  had 
a  turbulent  stack  gas. 

All  GLC  readings  were  obtained  with  a  moveable  sample  tube 
mounted  0.4  cm  (3.2  m  full  scale)  above  ground  level.  Four  consecutive 
samples  were  taken  at  each  downwind  position,  for  100  seconds  each. 

During  each  sample,  the  average  windspeed  was  measured  with 
a  Kurz  Model  435  anemometer  mounted  at  stack  height  beside  the  stack  at 
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TABLE  6-1:  SUMMARY  OF  GROUND  LEVEL  CONCENTRATION  MEASUREMENTS  WITH  PREDICTED  MAXIMUM  VALUES 
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From  predicted  values  by  (6-8) 
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x  =  0,  y  =  30  cm.  From  these  readings,  a  correction  was  made  to  each 
100  second  measurement  for  small  deviations  in  mean  flow  velocity  about 
the  desired  test  value  (see  Appendix  A).  This  correction  was  at  most 
about  10%,  for  cumulative  400  second  averages. 

Limitations  of  the  Study 

The  results  of  this  study  were  limited  by  two  factors: 

First,  the  tunnel  test  section  downwind  of  the  stack  was  not 
long  enough  to  locate  maximum  GLC  values  for  all  the  plumes  tested. 

This  was  particularly  true  of  the  two  low-speed  cases  for  plume  models 
A  and  B,  where  measured  concentrations  rose  to  only  about  half  of  the 
maximum  GLC  at  the  end  of  the  test  section,  as  shown  in  Figure  6-1. 

For  the  high  speed  cases  of  these  models,  the  position,  x  ,  of 
maximum  GLC  was  found  to  be  near  the  end  of  the  test  section. 

As  will  be  shown  later,  the  lack  of  far  downstream  data  only 

allowed  the  use  of  a  simple  model  for  dike  effects. 

The  second  limitation  of  the  study  involved  the  fact  that  the 
data  for  plume  model  A  was  inherently  less  useful  than  the  data  for 
plume  model  B.  This  was  because  during  all  model  A  runs,  a  cloth 
filter  was  installed  across  the  end  of  the  tunnel  test  section  to  help 
to  steady  the  tunnel  velocity.  Unfortunately,  the  filter  was  later 
found  to  alter  mean  flow  patterns  in  the  tunnel  so  that  the  plume 

gradually  drifted  toward  one  side  of  the  tunnel,  as  it  passed  through 

the  test  section. 

Because  GLC  measurements  were  taken  along  the  tunnel  center- 
line  (assuming  no  lateral  shift)  they  were  lower  than  the  values 
along  the  plume  axis.  This  was  investigated  by  finding  the  shifted 
maximum  GLC  line  by  measuring  lateral  profiles  and  taking  one  GLC 
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profile  down  this  line  for  each  of  the  two  velocities  with  flat  terrain. 
By  comparison  with  the  corresponding  flat  terrain  profiles  down  the 
tunnel  centerline,  there  was  no  appreciable  difference  at  Us  =  0.5  m/sec, 
but  a  significant  (20-50%)  difference,  increasing  with  x,  for  the 
profiles  at  =  1.0  m/sec.  The  measured  ratios  of  GLC  between  the 
maximum  and  tunnel  centerline  values  from  this  test  were  used  to 
correct  all  1  m/sec  model  A  data. 

The  result  of  the  above  was  that  while  model  B  data  could 
be  used  to  explore  all  aspects  of  absolute-value  prediction  of  GLC 
with  and  without  dike  effects,  the  model  A  data  was  uncertain  in  an 
absolute  sense  and  was  limited  only  to  the  examination  of  dike  effects 
relative  to  the  flat  terrain  case. 

Gaussian  Model  for  Plume  Dispersion 

The  Gaussian  plume  model  was  found  to  be  very  useful  in 
predicting  GLC  values.  Numerous  investigations,  such  as  the  full 
scale  studies  of  Weil  and  Jepsen  (1977)  and  laboratory  measurements 
by  Wilson  (1977),  have  shown  that  the  statistical  normal  or  Gaussian 
distribution  is  a  close  model  to  measured  concentration  distributions 
produced  by  turbulent  diffusion.  Applied  to  the  case  of  a  smokestack 
in  an  atmospheric  boundary  layer,  the  Gaussian  solution  for  plume 
mass  concentration  X  is: 


X(x,y,z,h) 


exp 


(6-1) 
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In  the  model,  m  is  the  mass  flow  rate  of  contaminant  from  the  stack  and 

a  and  a  are  the  standard  deviations  of  y  and  z  direction  concentra- 
y  z  ' 

tion  profiles,  the  "plume  spreads". 

Beneath  the  plume  centerline  at  ground  level  (z=y=0) ,  (6-1) 
reduces  to: 


m 


tto  a  U 
y  z 


(6-2) 


Using  m  =  XSQ  where  Q  is  the  total  stack  volume  flow  rate,  and  special¬ 
izing  to  the  case  of  isothermal,  isobaric  conditions  between  stack  and 
atmosphere,  such  as  exist  for  the  wind  tunnel  model,  (6-2)  may  be 


expressed  as  a  ratio  of  volume  concentrations: 


ttct  o  U 

y  z 


exp 


(6-3) 
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To  predict  GLC,  the  plume  height  h  and  the  plume  spreads  cr  and  a  must 

y  z 

be  specified  with  downwind  distance. 

In  this  study,  the  measured  average  value  of  windspeed  was 

used  with  (6-3).  Turner  (1970)  suggests  that  the  integrated  mean  value 

U  in  the  plume  layer  should  be  used,  which  is  not  much  different  than 

U  here, 
s 

Plume  Spread 

Lateral  and  vertical  plume  spreads  for  the  flat  terrain  cases 
of  plume  model  B  were  obtained  by  fitting  Gaussian  distributions  to 
measured  lateral  and  vertical  concentration  profiles.  For  the  vertical 
profiles,  the  fits  were  based  only  on  the  data  below  the  plume  center- 
line  at  h.  Vertical  plume  spreads  defined  in  this  way  should  be  more 
relevant  for  GLC  prediction  than  full-profile  fits.  An  example  of  a 
fitted  vertical  profile  is  shown  in  Figure  6-2.  It  can  be  seen  that 
the  lower  half  of  the  profile  was  well-fitted  by  a  (lower)  =  10.5  cm 

'Zj 

at  the  expense  of  a  slightly  poorer  fit  for  the  top  half.  The  fit  shown 
based  on  the  entire  profile,  with  o -  9.4  cm,  was  noticeably  different 
in  this  case,  although  for  most  profiles,  it  was  found  that  the  choice 
of  a fitted  to  lower-half  or  full  profiles  made  very  little  difference 
to  GLC  prediction. 

Because  it  was  necessary  to  interpolate  and  extrapolate  spread 
values  for  GLC  prediction,  a  smoothed  functional  form  was  required.  Full 
scale  dispersion  studies  such  as  the  Brookhaven  study,  Singer  and  Smith 
(1966),  have  shown  that  for  a  given  atmospheric  condition,  spreads  of 
a  non-buoyant  contaminant  in  a  given  direction  are  well  fitted  by  a 
single  power  law,  a  =  ax  .  To  reduce  measured  plume  spreads,  a  (total). 
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which  contain  the  effects  of  buoyancy-generated  spread  to  their  corres¬ 
ponding  values  in  non-buoyant  form,  a  (non-buoyant) ,  Winkel  (1979) 
suggests  the  use  of  the  equation 

2  2  2 
a  (total)  =  a  (non-buoyant)  +  a  (correction)  (6-4) 

where  the  latter  term,  representing  buoyancy- induced  spread,  may  be 
estimated  by  the  following: 

2  2  2 

O  (correction)  =  C  3  Ah  (6-5) 

In  this  correction  term,  and  Ah  are  those  values  which  would  apply 
to  the  rise  of  the  plume  in  a  laminar  cross-flow,  which  in  turn  is  re¬ 
lated  to  buoyancy- induced  spread.  The  constant  was  found  by  Winkel 
(1979)  to  be  about  1/9.  From  rise  measurements  of  plume  model  B  in  a 
laminar  cross-flow,  3 2  =  0.87,  and  this  value  was  used  for  buoyancy 
corrections  throughout  the  study. 

Using  (6-4)  and  (6-5),  the  measured  a  and  a  (lower)  values 

y  z 

for  plume  model  B  were  buoyancy  corrected  to  a (non-buoyant)  form.  In 
Figures  6-3  and  6-4,  power  law  fits  to  this  data  are  plotted  and  compared 
to  800:1  reduced  Brookhaven  C  (neutral  stability)  power  laws. 

To  obtain  plume  spread  values  for  any  downwind  position,  the 
fitted  power  laws  from  plume  model  B  were  corrected  back  to  model  condi¬ 
tions  using  (6-4)  and  (6-5).  Thus,  the  magnitude  of  the  buoyancy 
correction  was  not  very  important;  the  advantage  of  using  the  correction 
was  to  obtain  a  universal  functional  form  for  spread. 

Predicting  GLC  for  the  Flat  Terrain  Case 


Using  (6-3),  profiles  of  Cq/Cs  were  predicted  for  the  two  tested 
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windspeeds  for  plume  model  B,  as  plotted  in  Figure  6-5.  Rise  values 
used  in  the  calculations  were  obtained  from  (4-11)  to  match  the  actual 
measured  rise.  For  =  0.45m/s,  the  combined-rise  model  shown  in 
Figure  4-1  was  used,  with  32  =  1.4.  The  final  rise  was  set  at 
^final  =  16.37  cm  at  x£^na-^  =  2200Lg.  For  =  lm/sec,  Ah  =  0  was 
used  for  all  x  due  to  stack  downwash. 

Several  important  conclusions  may  be  drawn  from  Figure  6-5: 

(1)  The  difference  in  GLC  prediction  between  using  from  whole- 
profile  fits  or  a  (lower)  is  negligible.  Predicted  values 

Lj 

from  whole-profile  Gaussian  fits  are  very  close  to  those  using 
a  (lower)  =  0.271  x^'7(\ 

(2)  The  Gaussian  model  for  plume  dispersion  appears  to  be  valid 
under  a  wide  range  of  conditions,  from  the  case  with  severe 
downwash  to  those  with  large  amounts  of  plume  rise.  The 
vertical  and  lateral  measured  plume  profiles  used  to  obtain 
plume  rise  and  spread  were  well  represented  by  Gaussian  dis¬ 
tributions,  as  seen  in  Figure  6-2.  The  good  agreement  of  the 
measured  GLC  profiles  in  Figure  6-5  to  the  Gaussian  predictions 
also  supports  the  use  of  the  Gaussian  model  because  the  plume 
rise  and  spread  values  used  in  the  predictions  were  not  derived 
in  any  way  from  ground  level  data. 

Figure  6-5  also  displays  some  inadequacies  of  the  Gaussian 
model  for  GLC  prediction.  For  example,  there  is  a  kink  in  the  pre¬ 
dicted  curve  for  U  =0.45  m/s  near  x  =  225  cm,  and  the  curve  is  seen 

s 

to  rise  more  sharply  than  the  data  beyond  that  point.  The  scatter  in 
the  data,  which  was  typically  about  +0.01  for  (Cq/Cs)  x  1000,  can 
not  explain  this  discrepancy.  Actually,  the  change  in  character  of 
the  prediction  was  caused  by  the  final  rise  model,  where  xf ^nai=2200Lg 
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(223  cm) .  Final  rise  was  observed  for  this  plume,  as  seen  in  Figure 
4-1,  but  the  data  in  Figure  6-5  behaved  as  if  the  plume  did  not  level 
off.  Using  the  Gaussian  model,  plume  centerline  changes  were  assumed 
to  be  instantly  reflected  in  values  of  GLC,  however,  the  data  was 
accurate  enough  to  display  the  error  in  this  assumption.  There  is 
probably  a  time  delay  for  the  edges  of  a  plume  to  respond  to  centerline 
changes . 

Predicted  GLC  curves  for  the  plume  model  A  flat  terrain  cases 
were  also  required.  Because  no  rise  data  was  available,  and  because 
the  plume  model  A  data  was  already  limited  to  the  investigation  of 
relative  effects  of  dikes,  the  predictions  for  these  cases  were  fitted 
to  the  data,  as  shown  in  Figure  6-6,  by  the  choice  of  an  appropriate 
entrainment  constant,  32*  The  advantage  of  this  was  that  when  the  flat 
terrain  predictions  were  later  corrected  for  dike  effects  and  compared 
to  data,  the  closeness  of  fit  gave  a  direct  indication  of  the  accuracy 
of  the  dike  correction.  The  closeness  of  GLC  predictions  to  plume  model 
A  data  are  not  a  verification  of  the  absolute  accuracy  of  the  Gaussian 
model,  although  the  rise  values  used  to  obtain  the  fit  for  the  flat 
terrain  cases  were  very  realistic:  for  the  1  m/sec  case,  z  =  0  (all  x) 
was  used,  and  for  the  0.5  m/sec  case,  (4-11)  was  used  with  final  rise 
at  x  =  2200L-.  and  30  =  1.22  (compared  to  30  =  1.4  for  the  turbulent 
stack  gas  in  model  B) . 

Correcting  for  Dike  Effects 

A  model  to  correct  for  dike  effects  on  GLC  was  obtained  by 
taking  the  ratio  between  the  expressions  for  predicted  GLC,  with  and 
without  a  dike,  using  (6-3) : 
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Cq (dike) 
C 

o 


0^ (dike)  G^dike)  eX^ 


-h2 (dike) 

2  az2(dike) 


(6-6) 


In  (6-6),  the  parameters  C  ,  a  ,  a  ,  and  h  are  for  flat  ter- 

o  y  z 

rain.  Defining  some  "dike  disturbance"  parameters. 


n 


a 

z 


°z (dike) 


£ 


a 

y 

Gy (dike) 


h(dike) 

h 


(6-6)  becomes: 


Co (dike) 
C 

o 


ep  exp 


'  2  2 
Y  n 


(6-7) 


With  suitable  models  for  p,  e,  and  y,  (6-7)  would  account  for 
all  possible  dike  disturbance  effects  in  GLC  prediction  using  the  Gaus¬ 
sian  model.  However,  in  comparing  the  model  to  data,  it  was  practical 
to  allow  only  one  of  these  parameters  to  vary,  setting  the  others  equal 
to  unity  (assuming  no  effect  of  dike)  or  equal  to  the  varied  parameter 
(assuming  equal  dike  effect).  A  number  of  possibilities  were  examined, 
as  summarized  below: 


(1)  y  =  l,  e  =  n 

Assigning  dike  effects  equally  to  o  and  o  whil< 

y  z 


keeping  plume 
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rise  constant  seemed  physically  reasonable  but  did  not  produce  a  workable 
correction.  Equation  (6-7)  was  unable  to  predict  a  large  enough  cor¬ 
rection  effect  in  many  cases  regardless  of  the  q  values  chosen. 

(2)  e  =  1,  n  =  1 


Assigning  dike  effects  to  plume  rise  alone  while  keeping  plume 
spread  constant  always  produced  a  proper  correction  if  y  was  varied, 
because  (6-7)  is  unbounded  in  CQ(dike)/Co  as  y  decreases  from  y  =  1. 
However,  this  model  had  the  disadvantage  of  being  physically  unrealistic, 
considering  that  a  constant  height  non-buoyant  plume  which  would  not 
be  subject  to  correction  could  still  be  affected  by  a  dike. 

(3)  e  =  1,  y  =  1 

Assigning  dike  effects  to  a  alone  while  keeping  o  and  h 

z  y 

constant  was  found  to  be  both  workable  and  physically  reasonable.  Extra 
turbulence  induced  by  a  dike  wake  would  be  expected  to  increase  vertical 
spread  and  cause  larger  GLC  values,  which  was  consistent  with  the  model 
and  with  measurements.  Lateral  plume  spread  would  also  be  expected  to 
be  increased  by  dike-induced  turbulence,  which  causes  some  loss 
of  realism  in  the  model,  however,  this  did  not  cause  it  to  be  untenable 
as  in  (2)  above. 

With  the  use  of  £  =  1  and  y  =  1,  (6-7)  is: 


C  (dike) 
o 

C 

o 


(6-8) 


A  sketch  of  this  function  for  a  typical  value  of 


1.1 


is  shown  below: 
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Only  that  part  of  the  function  (shown  as  a  solid  line)  where  C  (dike)/C 

o  o 

increases  with  decreasing  p  was  physically  realistic  and  used  for  the 
correction;  the  other  part  (shown  as  a  dashed  line)  was  extraneous.  Note 
that  there  is  a  maximum  correction  possible  by  this  model  for  a  given 


h2/2az2 


value.  Measured  values  of  C  (dike)/C  were  never  found  to 

o  J  o 


exceed  this  maximum, 


Using  smoothed  curves  through  all  the  measured  GLC  profiles, 
values  of  CQ(dike)/Co  were  obtained  at  about  nine  regularly-spaced 
downwind  positions  for  each  case.  Then,  using  (6-8),  p  was  varied 
downwards  from  p  =  1  until  a  match  was  obtained  in  CQ(dike)/Co  for 
each  position.  The  final  values  of  p  required  to  match  experimental 
measurements  are  shown  as  data  points  in  Figure  6-7. 

It  was  expected  that  p  would  be  equal  to  unity  (no  dike  effect) 
upstream  of  the  dikes,  and  that  its  value  would  decrease  downwind  of 
the  dikes  and  then  return  to  unity  at  a  certain  number  of  dike  heights 
downwind  of  the  dikes  where  the  dike  wake  turbulence  had  decayed  to 
ambient  levels.  This  behaviour  was  observed  (with  a  great  deal  of  scatter) 
for  the  medium  (45.8  m  full  scale)  dike  data  in  Figure  6-7,  with  the 
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effect  of  the  dike  dying  away  at  about  80  dike  heights  downstream. 

For  the  large  (94  m  full  scale)  dike  the  limitation  of  test  section  length 
did  not  allow  this  return  to  r|  =  1  far  downstream  to  be  observed,  al¬ 
though  the  measured  ri  values  were  consistent  in  form  to  those  of  the 
medium  dike.  The  r|  values  for  the  large  dike  indicated  a  much  increased 
dike  disturbance  effect  over  the  medium  dike  values. 

The  low  values  of  p-0.7  observed  for  two  of  the  large  dike 
cases  in  Figure  6-7  occurred  as  Equation  6-8  was  approaching  a  maximum 
in  Co(dike)/CQ.  In  this  region,  the  correction  is  very  insensitive 
to  changes  in  p ,  and  ignoring  these  deviations  in  the  recommended 
line  for  p  shown  in  Figure  6-7  does  not  produce  large  errors  in  predicting 
the  dike  correction. 

A  Model  for  Vertical  Spread  Ratio,  p 

Owing  to  the  fairly  rough  and  incomplete  nature  of  the  data 
for  P,  only  a  simple  model  was  proposed  to  correlate  the  data  and 
allow  for  the  prediction  of  GLC  with  dikes.  Straight  line  segments 
were  used,  which  divided  the  model  into  five  regions  in  xd  =  [x-x(dike)J /H 
as  given  below  and  shown  in  Figure  6-7. 


x 


d 


<-2.5 


^low  -  1 
n  ”  Xd  22.5 


+ 


-2 . 5<  xd<  20 


20<  x,<  60 


(6-9) 
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p  =  xd 
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P  =  1 
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The  minimum  value  of  q  ,  which  defines  the  magnitude  of  the  correction, 
is  given  by 


n 


low 


1  - 


(6-10) 


s 

— 

It  is  clear  that  from  the  data  in  Figure  6-7  that  (1  -  q,  )  increases 
non-linearly  with  dike  height,  and  it  is  also  expected  that  the  ratio 
of  the  dike  height  to  the  stack  height  would  be  of  importance  in  a 
general  correlation.  Because  only  two  dike  heights  and  one  stack  height 
were  tested,  (6-10)  is  only  a  preliminary  estimate  in  lieu  of  further 
data. 


In  (6-9),  the  extension  of  dike  effects  slightly  upwind  of  the 
dike  crest  (2.5  dike  heights)  was  to  account  for  local  streamline  com¬ 
pression  effects  on  GLC  which  were  observed  to  increase  GLC  values 
directly  above  the  dikes.  This  phenomenon  was  also  observed  by  Wilson 
(1977)  for  flow  over  gentle  hills. 


Predicted  GLC  With  Dikes 

Equations  (6-8),  (6-9)  and  (6-10)  were  used  to  predict  GLC 
values  for  all  the  cases  tested  in  the  study  (see  Table  6-1).  These 
predictions  are  shown  along  with  the  measured  values  in  Figures  6-8 
through  6-12.  In  each  case,  the  predicted  curve  for  the  flat  terrain 
case  is  also  shown. 

A  number  of  conclusions  may  be  drawn  from  these  predictions: 

(1)  The  correction  model  was  successful,  producing  about  the  right 
amount  of  increase  in  GLC  above  the  flat  terrain  prediction  for  all 
cases  (considering  that  the  data  was  used  to  obtain  the  model  for  q  this 
is  not  too  surprising).  For  the  plume  model  A  cases,  this  can  be  directly 
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seen  from  the  good  closeness  of  fit  between  the  data  and  the  predicted 
curves  with  correction.  For  the  1  m/sec  case  of  plume  Model  B  shown 
in  Figure  6-12,  the  apparent  overprediction  in  the  dike  correction  was 
due  to  the  fact  that  the  flat  terrain  prediction  also  overpredicted  the 
flat  terrain  data  by  about  the  same  amount,  as  shown  in  Figure  6-5. 

In  Figure  6-12,  a  smoothed  curve  from  data  taken  with  all  condi¬ 
tions  the  same  except  the  plant  site  model  reversed  is  seen  to  follow 
the  data  for  =  0.9  m/sec  very  closely.  This  indicates  that  the  dike 
was  by  far  the  dominant  terrain  feature  and  that  results  from  this  study 
may  be  generalized  to  any  typical  plant  site. 

(2)  Dike  effects  on  GLC  are  definitely  not  linear  with  dike  height. 

The  large  dike  is  seen  in  Figures  6-8  through  6-12  to  produce  increases 
in  GLC  over  the  flat  terrain  case  about  an  order  of  magnitude  greater 
than  those  for  the  medium  dike,  which  was  about  half  the  size.  This 
is  accounted  for  in  the  model  in  Equation  6-10  where  (1  -  is 

3 

proportional  to  H  . 

Maximum  Values  of  Ground  Level  Concentrations  (GLC) 

For  all  the  cases  tested,  the  GLC  predictions  from  the  Gaussian 
plume  model  (6-8),  with  dike  corrections  (6-9)  and  (6-10)  were  carried 
out  as  far  downwind  as  necessary  to  determine  the  position,  xmax>  an<i 
magnitude  of  the  maximum  value  of  GLC.  These  results  are  given  in 
Table  6-1  where  the  magnitude  results  are  expressed  as  the  ratio  of 
maximum  GLC  with  the  dike  to  the  maximum  GLC  for  the  flat  terrain  case. 

In  Table  6-1  it  is  seen  that  the  effect  of  the  large  dike  was 
to  increase  the  maximum  GLC  value  by  about  20%  in  all  cases,  and  also 
cause  that  value  to  be  reached  about  20%  closer  to  the  stack  compared 
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to  the  flat  terrain  case.  The  medium  (45.8m  full  scale)  dike  was  found 
to  have  only  at  most  a  3%  effect  on  maximum  GLC,  with  a  negligible 
effect  on  x 

max 
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FIGURE  6-1:  PREDICTED  GLC  FOR  PLUME  MODEL  A  WITH  FLAT  TERRAIN 
AND  UPWIND  DIKES 


MODEL  HEIGHT  RBOVE  GROUND  (CM) 
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FIGURE  6-2:  GAUSSIAN  FITS  TO  A  VERTICAL  CONCENTRATION  PROFILE 
FROM  PLUME  MODEL  B  AT  x  =  150  cm 


z  (lower)  Model  Buoyancy-Corrected  Vertical  Plume  Spread  (cm) 
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FIGURE  6-3:  BUOYANCY-CORRECTED  VERTICAL  SPREADS  FROM  BEST-FIT 
TO  LOWER  HALF  OF  CONCENTRATION  PROFILES  FOR  PLUME 
MODEL  B 
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FIGURE  6-4:  BUOYANCY-CORRECTED  LATERAL  SPREADS  FOR  PLUME  MODEL  B 
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X,  Model  Scale  Downwind  Distance  from  Stack,  cm 
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X,  Full  Scale  Downwind  Distance  from  Stack,  meters 
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FIGURE  6-5 :  PREDICTED  AND  MEASURED  GLC  FOR  PLUME  MODEL  B  WITH 
FLAT  TERRAIN 
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X,  Model  Scale  Downwind  Distance  from  Stack,  cm 
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X,  Full  Scale  Downwind  Distance  from  Stack,  meters 
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FIGURE  6-6: 


PREDICTED  AND  MEASURED  GLC  FOR  PLUME  MODEL  A  WITH 
FLAT  TERRAIN 


Theory  (H  =  45.8  m) 
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FIGURE  6-7:  DIKE  DISTURBANCE  PARAMETER  n  FROM  MEASURED  GLC 
FOR  PLUME  MODELS  A  AND  B 
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X,  Model  Scale  Downwind  Distance  from  Stack,  cm 
800  1000  1200  1400  1600  1800  2000  2200  2400 

X,  Full  Scale  Downwind  Distance  from  Stack,  meters 
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FIGURE  6-8:  PREDICTED  AND  MEASURED  GLC  FOR  PLUME  MODEL  A  AT 

U  =  1  0  m/s  WITH  THE  MEDIUM  DIKE  UPWIND  AND  DOWNWIND 
s 
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X,  Model  Scale  Downwind  Distance  from  Stack,  cm 
800  1000  1200  1400  1600  1800  2000  2200  2400 

X,  Full  Scale  Downwind  Distance  from  Stack,  meters 
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X,  Model  Scale  Downwind  Distance  from  Stack,  cm 
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X,  Full  Scale  Downwind  Distance  from  Stack,  meters 
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FIGURE  6-9:  PREDICTED  AND  MEASURED  GLC  FOR  PLUME  MODEL  A  AT 

U  =  1.0  m/s  WITH  THE  LARGE  DIKE  UPWIND  AND  DOWNWIND 
s 
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X,  Model  Scale  Downwind  Distance  from  Stack,  cm 

800  1000  1200  1400  1600  1800  2000  2200  2400 

X,  Full  Scale  Downwind  Distance  from  Stack,  meters 
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FIGURE  6-10: 


PREDICTED  AND  MEASURED  GLC  FOR  PLUME  MODEL  A  AT 
U  =0.5  m/s  WITH  THE  MEDIUM  DIKE  UPWIND  AND  DOWNWIND 
s 
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X,  Full  Scale  Downwind  Distance  from  Stack,  meters 
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FIGURE  6-11:  PREDICTED  AND  MEASURED  GLC  FOR  PLUME  MODEL  A  AT 

U  =0.5  m/s  WITH  THE  LARGE  DIKE  UPWIND  AND  DOWNWIND 
s 
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FIGURE  6-12:  PREDICTED  AND  MEASURED  GLC  FOR  PLUME  MODEL  B  WITH 
DIKES  UPWIND  AND  DOWNWIND 


CHAPTER  VII 


SUMMARY  AND  RECOMMENDATIONS 


A  wind  tunnel  model  was  used  to  examine  the  behaviour  of 
stack  plumes  on  an  oil  sands  surface-mining  site.  As  a  result  of 
the  current  modelling  practises,  the  mean  flow  Reynolds  number  in 
the  wind  tunnel  was  about  25,000  times  smaller  than  that  in  the 
full  scale  atmospheric  flow.  This  required  a  careful  investigation 
of  Reynolds  number  effect  on  several  flow  situations  in  the  model 
before  the  model  could  be  used  with  confidence.  The  results  of  these 
investigations  are  applicable  to  many  physical  modelling  situations 
where  Reynolds  number  mismatches  occur. 

Reynolds  Number  Effects  on  Cylinder  Flows 

The  wakes  behind  the  model  stacks  appear  to  exert  realistic 
downeffects  on  the  model  plumes  in  spite  of  Reynolds  number  mismatch. 
This  was  due  to  the  close  similarity  in  cylinder  base  pressure  coeffi¬ 
cients  at  -C„.  -0.86  between  the  model  stack  flows  in  the  subcritical 
Pb 

regime  and  the  full  scale  stack  flows  in  the  transcritical  regime. 

Great  care  must  be  taken  in  modelling  of  the  stack  wake  effect,  because 
if  the  full  scale  stack  flows  are  in  the  supercritical  regime,  where 
base  pressure  coefficients  -Cpb  are  much  smaller,  the  model  stacks  will 
likely  exaggerate  stack  downwash. 

Plume  Rise  and  Stack  Downwash 

A  plume  rise  model  for  neutrally  stable  atmospheres,  (4-11), 
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was  developed  and  included  the  effects  of  both  plume  buoyancy  and 
momentum.  This  model  was  found  to  compare  very  well  with  measured 
plume  trajectories.  It  appears  that  momentum  rise  (the  1/3  law)  and 
buoyancy  rise  (the  2/3  law)  add  as  the  sum  of  cubes.  Therefore,  either 
one  is  equally  as  useful  a  means  of  producing  plume  rise  as  the  other. 
Also,  for  the  majority  of  the  plume  trajectory,  one  type  of  rise  is 
dominant,  making  it  reasonable  to  consider  momentum  jet  entrainment 
and  buoyant  plume  entrainment  as  separate  phenomena. 

Values  of  the  momentum  rise  entrainment  constant  3^  to  match 
measured  plume  rise  were  predictable  from  a  simple  semi-empirical 
model,  (4-14).  This  model  is  recommended  for  use  in  general  stack 
design  calculations. 

The  buoyancy  rise  entrainment  constant  39  varied  with  atmos¬ 
pheric  and  stack  flow  conditions  in  agreement  with  physical  arguments 
regarding  entrainment.  It  was  possible  to  identify  constant  3 2 
values  for  various  flow  conditions,  although  the  validity  of  such  a 
procedure  is  somewhat  questionable.  Apparently,  there  is  no  univer¬ 
sally  applicable  constant  value  of  3  ,  as  suggested  by  Briggs  (1975) 
(3  =0.6).  Developing  a  comprehensive  model  for  39  to  include 

all  relevant  entrainment  effects  is  a  topic  for  future  study. 

Plumes  in  turbulent  cross-flows  were  observed  to  achieve  final 
rise.  The  simple  criteria,  xfinal  =  2200  LR,  was  found  to  predict 
this  phenomenon  for  plumes  under  a  wide  variety  of  conditions.  The 
constant  of  proportionality  might  have  to  be  adjusted  slightly  for  use 
with  boundary  layers  other  than  the  one  tested. 

The  momentum  flux  ratio  4^  was  found  to  be  strongly  correlated 
with  the  degree  of  severity  of  the  stack  downwash  effect.  To  avoid 
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downwash,  the  designer  should  ensure  that  d>w  >  2.  At  lower  values, 

M 

plume  rise  will  be  suppressed,  but  still  predictable  by  (4-11)  down  to 
about  =  1,  using  adjusted  S2  values.  Below  <j>  =  1,  little  or 

no  plume  rise  is  expected. 

Reynolds  Number  Effects  on  Dike  Flows 

Flows  past  tailings  pond  dikes  were  found  to  be  severely  affected 
by  the  Reynolds  number  mismatch  in  the  model.  Model  dikes  had  a 
greatly  exaggerated  wake  size  and  influence  on  downwind  flow  turbulence. 
This  type  of  modelling  error  is  expected  for  all  types  of  terrain 
features  where  separation  points  are  not  fixed  by  sharp  edges.  The 
use  of  a  flow  deflector  vane  mounted  on  the  dike  crest  improved  the 
model  substantially,  creating  mean  flow  and  turbulence  profiles  down¬ 
wind  of  the  dikes  which  were  close  to  the  expected  full  scale  profiles. 

Dike  Effects  on  Ground  Level  Concentration  (GLC) 

The  Gaussian  model  for  plume  dispersion  was  found  to  be 
adequate  in  fitting  both  measured  plume  concentration  profiles  and 
predicting  GLC  values.  More  complex  dispersion  models  do  not  seem  to 
be  necessary,  as  the  Gaussian  model  could  predict  all  the  important 
plume  behaviour,  such  as  ground  reflection,  and  the  shape  of  ground 
level  concentration  profiles. 

The  presence  of  dikes  upwind  and  downwind  of  a  stack  was  found 
to  increase  measured  GLC  values  in  general.  This  could  be  accounted 
for  in  the  Gaussian  model  by  the  use  of  increased  a  values.  In- 

Lt 

creases  in  a  are  not  the  only  dike  disturbance  effect,  but  the  fact 
z 

that  this  was  the  one  model,  of  the  attempted  corrections  that  gave 
reasonable  results  suggests  that  the  dominant  effect  of  dikes  may  be 
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on  a  . 
z 

Dike  effects  on  GLC  were  found  to  be  very  non-linear  with  dike 

3 

height,  as  provided  for  by  the  dependence  of  q  on  H  in  (6-10),  and 
the  non-linearity  of  the  Gaussian  model  itself,  in  (6-8).  Thus, 
large  dikes  produced  disproportionately  larger  disturbances  to  GLC 
than  smaller  dikes. 

The  maximum  value  of  GLC  for  a  plume  was  found  to  be  fairly 
insensitive  to  dike  disturbances,  even  when  large  increases  in  GLC 
due  to  dikes  were  present  prior  to  x  .  For  dikes  of  about  half 
the  stack  height  in  the  present  model,  increases  in  C0max  were  only 
about  20%.  This  means  that  predicting  C0max  assuming  flat  terrain 
would  in  most  cases  be  adequate  for  stack  design. 

Several  specific  recommendations  are  made  from  the  results 
of  the  present  study: 

(1)  In  plume  modelling,  models  should  be  as  large  as  conveniently 
possible  to  minimize  Reynolds  number  mismatch. 

(2)  When  modelling  stack  downwash,  smooth  model  stacks  may  probably 
be  used  if  the  full  scale  stack  flows  are  transcrit ical .  However, 
full-scale  base  pressure  measurements  on  stacks  would  be  very  helpful 
in  determining  if  the  flows  were  correctly  modelled,  especially 
because  of  the  stack  end  effect.  Further  investigation  in  this  area 
is  recommended. 

(3)  Stack  designers  should  consider  momentum  as  well  as  buoyancy 
as  a  means  of  producing  plume  rise,  particularly  because  increasing 

d)w  will  minimize  stack  downwash. 

M 

(4)  When  modelling  flow  over  smooth  terrain  obstacles  such  as  dikes, 
external  flow  modification  is  recommended  as  a  means  of  avoiding  the 
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unrealistically  large  recirculation  zones  that  typically  exist  at 
low  Reynolds  number. 

(5)  The  Gaussian  plume  model  is  recommended  for  use  in  predicting 
ground  level  concentrations  from  stacks.  The  far  downstream  variation 
in  ground  level  concentrations  should  be  investigated  to  determine 
if  the  Gaussian  model  will  remain  valid  after  a  buoyant  plume  has 
achieved  final  rise. 


119 


REFERENCES 


REFERENCES 


Arie,  M.  and  H.  Rouse  (1956),  "Experiments  on  Two  Dimensional  Flow 
Over  a  Normal  Wall",  J.  Fluid  Mech.  _1,  part  2,  pp  129-141. 

Briggs,  Gary  A.  (1969),  Plume  Rise,  USAEC  Critical  Review  Series, 

TID  -  25075,  Clearinghouse  for  Federal  Scientific  and  Technical 
Information. 

Briggs,  Gary  A.  (1970),  "Some  Recent  Analyses  of  Plume  Rise  Observa¬ 
tions",  Int.  Air  Poll.  Conference,  Wash.  D.C.,  Dec.  1970. 

Briggs,  Gary  A.  (1972),  "Discussion  of  Chimney  Plumes  in  Neutral  and 
Stable  Surroundings",  Atmos.  Environ.  6^,  pp  507-510. 

Briggs,  Gary  A.  (1975),  "Plume  Rise  Predictions",  Environmental 

Research  Laboratories,  Air  Resources,  Atmospheric  Turbulence 
and  Diffusion  Laboratory,  Oak  Ridge,  Tennessee,  June  1975, 

ATDL  Contribution  File  No.  75/15. 

Bringfelt,  B.  (1969),  "A  Study  of  Buoyant  Chimney  Plumes  in  Neutral 
and  Stable  Atmospheres",  Atmos.  Environ.  3_,  pp  609-623,  plus 
author's  reply  to  discussion  of  plume  rise  measurements  at 
industrial  chimneys,  pp  317-319. 

Bursnall,  W.J.  and  L.K.  Loftin  (1951),  "Experimental  Investigation 
of  Localized  Regions  of  Laminar  Boundary  Layer  Separation", 

Nat.  Adv.  Comm.  Aero.,  Wash.,  Technical  Note  2338. 

Chang,  S.C.  (1966),  "Velocity  Distributions  in  the  Separated  Flow 

Behind  a  Wedge  Shaped  Model  Hill",  Colorado  State  University 
Fluid  Dynamics  and  Diffusion  Laboratory  Technical  Report  CER 
65SCC66 . 

Counihan,  J.  (1969),  "An  Improved  Method  of  Simulating  an  Atmospheric 
Boundary  Layer  in  a  Wind  Tunnel",  Atmos.  Environ.  _3,  pp  197-214 

Counihan,  J.  (1975),  "Adiabatic  Atmospheric  Boundary  Layers:  A 
Review  and  Analysis  of  Data  from  the  Period  1880-1972", 

Atmos.  Environ.  9,  pp  871-905. 

Csanady,  G.T.  (1961),  "Some  Observations  on  Smoke  Plumes",  Int.  Jour. 
Air  Water  Pollut.  4_,  pp  47-51. 

Eliseev,  V.S.  (1973),  "Stereophotogrammetric  Investigation  of  the 

Air  Flow  in  the  Boundary  Layer  of  the  Atmosphere  Above  a  Hill", 
from  Berlyand,  M.  (ed.)  Air  Pollution  and  Atmospheric  Diffusion 
John  Wiley  §  Sons,  pp  95-108. 

Fage,  A.  and  V.M.  Falkner  (1931),  "The  Flow  Around  a  Circular  Cylinder' 
Aero.  Res.  Counc.,  London,  Rep.  and  Mem.  no.  1369. 


121 


Fay,  J.A.,  Escudier,  M.  and  D.P.  Hoult  (1969),  "A  Correlation  of  Field 
Observations  of  Plume  Rise",  Jour.  Air  Poll.  Control  Assn.  20, 
pp  391-397,  1970. 

Flachsbart,  0.  (1929),  from  an  article  by  H.  Muttray,  1932,  Handbook 
Experimental  -  Physik  £,  part  2  (Leipzig),  316. 

Good,  M.C.  and  P.N.  Joubert  (1968),  "The  Form  Drag  of  Two  Dimensional 
Bluff  Plates  Immersed  in  Turbulent  Boundary  Layers",  Jour. 

Fluid  Mech.  31 ,  part  3,  pp  547-582. 

Griffin,  Owen  M.  (1977),  "A  Universal  Strouhal  Number  for  the  'locking- 
on'  of  Vortex  Shedding  to  the  Vibrations  of  Bluff  Cylinders", 

Jour.  Fluid  Mech.  85,  part  3,  pp  591-606. 

Hinze,  J.O.  (1975),  Turbulence,  2nd  Edition,  McGraw  Hill,  p  643,  p  729. 

Isumov,  N.,  Jandali,  T.,  and  A.G.  Davenport  (1976),  "Model  Studies 

and  the  Prediction  of  Full  Scale  Levels  of  Stack  Gas  Concentration", 
Jour.  Air  Poll.  Control  Assn  26_,  pp  956-964. 

Ludwig,  G.R.  and  G.T.  Skinner  (1976),  "Wind  Tunnel  Modelling  Study 
of  the  Dispersion  of  Sulphur  Dioxide  in  Southern  Allegheny 
County,  Pennsylvania",  U.S.  EPA  Report  903/9-75-019,  pp  3-16. 

Moussa,  Z.M.,  Trischka,  John  W.  and  S.  Eskinazi  (1977),  "The  Near- 
Field  in  the  Mixing  of  a  Round  Jet  with  a  Cross-Stream", 

Jour.  Fluid  Mech.  80,  part  1,  pp  49-80. 

Patrick,  M.A.  (1967),  "Experimental  Investigation  of  the  Mixing  and 
Penetration  of  a  Round  Turbulent  Jet  Injected  Perpendicularly 
into  a  Transverse  Stream",  Trans.  Inst.  Chem.  Engrs.  (London) 

45,  pp  16-31. 

Plate,  E.J.  and  C.M.  Lin  (1964),  "The  Velocity  Field  Downstream  from 

a  Two-Dimensional  Model  Hill",  Part  1,  Colorado  State  University, 
Fort  Collins, Colo. ,  Technical  Report  No.  CER65EJP-CWL41 . 

Ricou,  F.P.  and  D.B.  Spalding  (1961),  "Measurements  of  Entrainments 
by  Axi symmetrical  Turbulent  Jets",  Jour.  Fluid  Mech.  £1, 
pp  21-32. 

Roshko,  Anatol  (1960),  "Experiments  on  the  Flow  Past  a  Circular 
Cylinder  at  Very  High  Reynolds  Number",  Jour.  Fluid  Mech. 

10,  pp  345-356. 

Sachs,  Peter  (1972),  Wind  Forces  in  Engineering,  2nd  Edition,  Pergamon 
Press,  1978,  p  70. 

Singer,  Irving  A.  and  Maynard  E.  Smith  (1966),  "Atmospheric  Dispersion 
’at  Brookhaven  National  Laboratory",  Int.  Jour,  of  Air  and  Water 
Pollution  10,  Pergamon  Press,  pp  125-135. 


122 


Slawson,  P.R.  (1978),  "Observations  and  Predictions  of  Natural 

Draft  Cooling  Tower  Plumes  at  Paradise  Stream  Plant",  Atmos. 
Environ.  12,  pp  1713-1724. 

Taylor,  G.I.  (1945),  "Dynamics  of  a  Mass  of  Hot  Gas  Rising  in  the 

Air",  USAEC  Report  MDDC-919  (LADC-276),  Los  Alamos  Scientific 
Lab. 


Weil,  Jeffrey  C.  and  Anders  F.  Jepsen  (1977),  "Evaluations  of  the 
Gaussian  Plume  Model  at  the  Dickerson  Power  Plant",  Atmos. 
Environ.  n_,  pp  901-910. 

Wilson,  D.J.  (1977),  "Plume  Diffusion  over  Downwind  Two-Dimensional 
Hills:  A  Wind  Tunnel  Study",  final  report  of  Alberta  Environ¬ 
ment  Research  Project  #75-8,  Department  of  Mechanical 
Engineering,  University  of  Alberta,  Edmonton,  Alberta. 

Wilson,  D.J.,  Winkel,  G.,  and  0.  Neiman  (1979),  "Reynolds  Number 

Effects  on  Flow  Recirculation  Behind  Two-Dimensional  Obstacles 
in  a  Turbulent  Boundary  Layer",  Proceedings  of  Fifth  Interna¬ 
tional  Wind  Engineering  Conference,  July,  1979,  Fort  Collins, 
Colorado. 

Winkel,  Gordon  (1979),  M.Sc.  Thesis,  Department  of  Mechanical 
Engineering,  University  of  Alberta,  Edmonton,  Alberta. 


APPENDIX  A 


124 


APPENDIX  A 

CORRECTION  TO  GLC  MEASUREMENTS  FOR  MEAN 
TUNNEL  SPEED  VARIATIONS 


In  the  present  study,  the  mean  tunnel  speed  was  found  to  vary 
slowly  around  the  desired  constant  value  for  a  test  by  about  +2% 

(maximum)  for  successive  100  second  averages.  It  was  possible  to  correct 
for  this  in  the  case  of  GLC  measurements. 

Ground  level  concentration  readings  are  very  sensitive  to 
changes  in  windspeed.  This  can  be  seen  by  examining  the  Gaussian  model 
for  GLC 


Q 


ttUo  a 
y  z 


(6-3) 


By  the  2/3  law  for  buoyant  plume  rise,  Ah  is  inversely  proportional 
to  U.  With  this,  (6-3)  becomes 


This  strong  dependence  of  Cq  on  windspeed  indicates  the  need  to  closely 
control  tunnel  speed  when  taking  GLC  measurements.  This  led  to  the  use 
of  the  cloth  filter  stretched  across  the  tunnel  test  section  end  for 
plume  model  A  measurements,  as  discussed  in  Chapter  VI. 

It  is  desirable  to  relate  small  changes  in  U  to  their  resulting 
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effect  on  C  .  Defining  a  concentration  parameter,  ip  = 
(6-3)  becomes 


(7rayazCo)/(QCs) , 


(A— 1 ) 


Taking  the  logarithm  of  (A-l), 


lnip  = 


(A-2) 


Using  small  disturbance  theory,  variations  are  assumed  only  in  ip  and 
U.  With  Ah  1/U,  (A-2)  is  differentiated  to  obtain 


dtjj  _  dU 
ip  ~  ~  U 


(Ah)h 
a  2 
z 


(A-3) 


Without  the  term  involving  Ah,  (A-3)  represents  the  relationship 

between  ip  and  U  for  a  non-buoyant  plume.  For  non-buoyant  sources, 

the  quantity  CU/Q  is  preserved,  and  this  differentiates  to  d\pAp  =  -  dU/U. 

2 

The  additional  term,  (Ah)h/u  ,  represents  the  modification  to  the 

z 

non-buoyant  relationship  to  account  for  plume  rise. 

Assuming  the  use  of  the  Brookhaven  C  (neutral  stability)  power 

0  78 

law  equation  for  vertical  spread,  0_  x  .  Also,  at  a  fixed  location, 

Z  " 

dip/ip  =  dCQ/Co  and  (A-3)  may  be  rewritten  as 


dU 

U 


K(Ah)h 

1.56 

x 


(A-4) 


where  K  is  a  dimensional  constant.  It  was  necessary  to  measure  K 
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because  the  amount  of  variation  measured  in  Cq  and  U  were  dependent 
on  the  way  measurements  were  taken,  through  sampling  times  used  and 
response  characteristics  of  the  equipment. 

A  study  was  undertaken  using  plume  model  A  at  Us  =  0.5  ^  0.1  m/s. 
Ground  level  concentration  measurements  were  taken  at  a  fixed  location, 
x  =  260  cm,  using  100  second  samples.  All  procedures  and  equipment 
were  as  described  in  Chapter  II  for  GLC  measurement. 

A  measurement  of  mean  windspeed  was  taken  before  and  after 
each  100  second  sample  using  the  Kurz  Model  435  Anemometer  with  capa¬ 
citor  damped  output  (time  constant  T  =  100  sec.).  The  average  of  these 
two  readings  was  used  as  an  estimate  of  the  average  value  of 
during  each  sample  period.  This  is  the  same  procedure  used  during  all 
GLC  measurement  in  the  present  study. 

From  a  large  number  of  readings  at  x  =  260  cm  at  windspeeds 

between  about  U  =0.4  m/s  and  U  =0.6  m/s  the  constant  K  was  ob- 

s  s 

tained  as  K  =  51  cm-0'44.  This  value  was  used  with  (A-4)  to  correct 
all  GLC  measurements  in  the  present  study.  For  example,  for  a  GLC 
sample  nominally  at  Us  =  0.5  m/sec,  the  correction  equation  would  be: 

Cq (measured)  -  Cq (standard)  0.5  -  Unmeasured) 

C  (standard)  0.5 

o 

where  C  (standard),  the  corrected  value  of  GLC,  is  the  only  unknown. 

Plume  rise  values  used  in  the  calculation  of  K  and  in  all 

model  A  corrections  were:  for  U  =  1  m/s,  Ah  =  2  cm,  and  for  U  =  0.5  m/s, 

s  ^ 

as  given  below: 


x  (cm) 

100 

110 

120 

130 

140 

150 

160 

170 

180  (and  up) 

Ah (cm) 

9.9 

10.4 

11.0 

11.6 

12.1 

12.5 

12.9 

13.4 

13.6 

K(Ah)h 

1.56 

x 
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When  (A-4)  was  used  to  correct  plume  Model  B  data,  Ah  =  0 

was  used  for  U  =  0.9  m/s  and  for  U  =0.45  m/s.  Ah  was  calculated 

s  s 

using  (4-11)  and  (4-15). 

The  magnitude  of  the  GLC  correction  using  (A-4)  was  typically 

less  than  10%  for  =  0.45  and  0.50  m/s,  and  typically  less  than 

1%  for  U  =0.90  and  1.0  m/s. 
s 


